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Preface 


The  purpose  of  this  study  was  to  compare  several 
techniques  for  the  prediction  of  mass-loaded  natural 
frequencies  and  mode  shapes.  Of  special  interest  in 
this  study  was  the  recovery  of  the  unloaded  mass,  stiff¬ 
ness,  and  damping  matrices  from  measured  modal  data 
using  a  non-square  modal  matrix  and  subsequent  solution 
for  mass- loaded  modal  data. 

This  study  was  somewhat • limited  in  scope  in  that 
only  one  panel  with  three  discrete  mass  loadings  was 
experimentally  tested  and  analysed.  It  is  hoped  the 
technique  using  pseudoinverses  will  be  explored  further 
to  determine  the  general  validity  of  this  method.  Appen¬ 
dices  A,  B,  and  C  should  be  helpful  in  this  endeavour. 

I  would  like  to  thank  my  advisor,  Capt.  H.  C,  Briggs 
of  the  Air  Force  Institute  of  Technology,  Mr.  R.  D.  Talmadge 
of  the  Air  Force  Flight  Dynamics  Laboratory,  and  Dr, 

P.  W.  Whaley  of  the  University  of  Nebraska  for  their 
support  and  guidance  in  this  effort.  Additionally,  I 
would  like  to  thank  my  wife,  Janice,  for  her  constant  en¬ 
couragement  and  inspiration.  Finally,  and  most  of  all, 

I  would  like  to  thank  my  Savior,  Jesus  Christ,  lor  that 
"Peace  which  passeth  all  understanding  (Phillipians  4:7)." 

In  keeping  with  Proverbs  3,  versus  5  and  6, 
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Abstract 


The  purpose  of  this  investigation  was  to  compare  the 
results  obtained  from  three  modal  prediction  techniques. 
The  first  technique  was  an  algorithm  developed  by  Whaley 
for  lightly  damped  structures  (Method  1).  Results  using 
this  algorithm  were  extracted  from  a  thesis  by  Glenesk. 

The  second  method  was  the  finite  element  method  using 
NASTRAN  (Method  2).  The  final  method  was  the  recovery  of 
unloaded  mass  and  stiffness  matrices  from  the  general  ma¬ 
trix-vector  differential  equation  of  modal  analysis  using 
modal  data  obtained  from  an  unloaded  test  item  (Method  3). 
Once  these  matrices  had  been  recovered,  a  quantity  of  mass 
was  added  to  the  mass  matrix  to  simulate  a  mass-loaded 
case.  The  generalized  eigenvalue  problem  was  solved  for 
mass-loaded  frequencies  and  mode  shapes  which  were  com¬ 
pared  to  experimental  results  for  the  same  test  item. 

Both  square  and  rectangular  modal  matrices  were  consider¬ 
ed  in  Method  3.  The  same  test  item  and  three  discrete 
mass-loaded  configurations  which  Glenesk  used  were  tested. 
Percentage  frequency  deviations  from  the  unloaded  test 
item  to  the  mass- loaded  predictions  ranged  from  -7.2-°d 
to  +7.4-%  in  Method  1,  from  -20.3-?4  to  +17,84-%  in  Method 
2,  and  from  -20.6-%  to  +8.4-%  in  Method  3.  Several  dis¬ 
crepancies  in  each  technique  prevent  a  direct  comparison 
of  these  results.  The  most  noteworthy  discrepancy  was 
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the  fact  that  the  modal  measurement  procedure  generated 
nonorthogonal  modes.  The  first  method  assumed  the 
mode  shapes  to  be  unaltered  between  the  unloaded  and  mass- 
loaded  cases  while  the  second  method  generated  mutually 
orthogonal  modes .  The  unloaded  nonorthogonal  mode  vectors 
were  used  in  Method  3  to  generated  mass-loaded  modal  quan¬ 
tities.  Detailed  procedures,  results,  and  conclusions 
are  obtained  in  the  body  and  appendices  of  the  report. 
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A  COMPARISON  OF  V.^RIOUS  TECHNIQ'/:.^  FOR 
THE  PREDICTION  OF  MASS-LOADED  MODE 
SHAPES  AND  NATURAL  FREQUENCIES 


I  Introduction 


Background 

The  ever- expanding  performance  envelopes  of  today's 
highly  complex  fighter  aircraft  subjects  them  to  increas¬ 
ingly  severe  vibration  environments.  Coupled  with  these 
severe  vibration  environments  is  a  desire  to  rapidly  in¬ 
corporate  newly-developed  weapons  system  technologies  into 
the  existing  fleet  of  fighter  aircraft.  One  such  example 
is  the  application  of  laser  physics  technology  to  vibration- 
sensitive  electro-optical  equipment  which  would  subsequent¬ 
ly  be  installed  in  high-performance  fighter  aircraft. 

The  installation  of  electro-optical  equipment  in  an 
aircraft  presents  a  complex  design  problem  in  that  it  is 
necessary  to  know  the  post-installation  modes  of  vibration 
and  natural  frequencies  of  the  aircraft  prior  to  the  actual 
installation  of  this  hardware.  Frequently,  the  only  modal 
data  available  to  the  designer  are  the  pre-installation 
modal  data.  Thus,  due  to  the  vibration  sensitivity  of 
the  electro-optical  hardware,  the  designer  must  consider 
how  to  properly  utiliiie  the  pre-installation  vibration 
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data  to  corroctly  predict  the  post- in.  i  Jlation  mode.s  o; 
vibration,  damping,  and  natural  frequencies,  Whaley  (Ref 
12)  summarized  three  analytical  techniques  to  accomplish 
this  task.  Additionally,  Glenesk  (Ref  3)  utilized  an 
algorithm  developed  by  Whaley  for  lightly  damped  struc¬ 
tures  (Ref  13)  to  predict  the  influence  of  added  lumped 
masses  on  the  vibration  characteristics  of  unloaded  struc¬ 
tures  . 

According  to  the  research  conducted  by  Glenesk  (Ref 
3),  as  the  size  of  the  added  mass  increased  relative  to 
the  mass  of  the  unloaded  structure,  the  accuracy  of  Whaley's 
algorithm  was  significantly  affected.  This  degradation 
in  algorithm  performance  might  be  attributed  to  the  assump¬ 
tion  that  the  unloaded  mode  shapes  are  unaffected  by  the 
addition  of  the  lumped  mass.  Another  possibility  involves 
the  fact  that  the  effect  of  damping  was  ignored  in  this 
algorithm.  Thus,  an  added  mass  might  have  significantly 
contributed  to  the  overall  structural  characteristics  in 
such  a  way  as  to  modify  the  mode  shapes  and  natural  fre¬ 
quencies  . 

In  addition  to  the  various  analytical  techniques,  a 
numerical  technique,  finite  element  analysis,  has  been 
widely  used  to  predict  mass-loaded  natural  frequencies, 
damping  ratios,  and  mode  shapes.  This  technique  requires 
the  construction  of  a  computer  model  in  which  the  continuous 


structure  is  idealized  as  a  conibinatio'  a  finite 


number  of  various  structural  components  (i.e.,  beams, 
rods,  plates,  etc.).  Although  accurate  results  can  be 
obtained  using  the  finite  element  method,  one  problem 
with  this  technique  is  the  significant  expenditure  of 
human  and  computer  resources  necessary  to  build,  debug, 
and  run  the  finite  element  code. 

The  rapid  development  of  portable  modal  analysis 
equipment  has  made  it  possible  to  eliminate  the  con¬ 
struction  of  a  finite  element  model  entirely.  With 
this  equipment  one  can  lay  out  a  suitable  grid  on  the 
portion  of  structure  to  be  modified,  conduct  standard 
modal  analysis  tests,  and  reduce  the  data  so  obtained 
to  determine  the  desired  unloaded  modal  data.  The 
question  then  becomes  how  to  properly  use  this  data 
to  determine  the  mass-loaded  modal  quantities  for  the 
modified  structure.  One  approach  to  this  dilemma  has 
been  suggested  by  Briggs  and  Whaley  (Ref  1)  whereby  one 
uses  the  general  matrix-vector  differential  equation  of 
structural  analysis,  its  solution  using  a  generalized 
coordinates  approach,  and  the  resulting  definitions 
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to  analytically  determine  the  mass-lc  ,  >1  modal  data 

when  only  the  experimentally  determined  unloaded  modal 
data  are  known. 

At  the  present  time,  the  usual  solution  to  equa¬ 
tions  1  through  3  requires  that  the  modal  matrices,  [u] 
T 

and  [U]  ,  be  square  matrices.  If  the  mass,  damping, 
and  stiffness  matrices,  [m],  [C],  and  [k],  are  n  x  n 
matrices,  [u]  and  [u]  must  also  be  n  x  n  matrices 
(where  "n"  is  the  number  of  measurement  points).  If 
one  measures  fewer  than  "n"  modes  in  the  frequency 
range  of  interest,  he  must  either  extend  this  frequen¬ 
cy  range  to  accomodate  "n”  modes,  or  reduce  the  grid 
size  to  "n"  grid  points.  As  either  of  these  approaches 
may  be  undesirable,  a  third  approach  using  the  method 
of  pseudoinverses  introduced  by  Penrose  (Ref  10)  may  be 
used  to  isolate  the  mass,  stiffness,  and  damping  ma¬ 
trices  on  the  left  hand  side  of  equations  1  through  3, 
respectively.  Appendix  A  contains  a  sample  problem 
for  the  reader  who  is  unfamiliar  with  this  technique. 
The  resulting  solution  will  be  an  approximate  solution 
to  the  mass,  damping,  and  stiffness  matrices  for  the 
unloaded  structure.  Then,  to  find  the  corresponding 
matrices  for  the  ma^s-loaded  structure  one  would  add 
appropriate  mass,  damping,  and  stiffness  quantities  at 
the  proper  locations  in  their  respective  matrices  to 
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simulate  the  structural  modification,  ird  resolv^e  equa¬ 
tions  1  through  3  for  the  mass-loaded  modal  information. 
If  damping  is  not  a  factor  one  wishes  to  consider  in 
this  analysis,  one  need  only  consider  the  solution  to 
the  standard  eigenvalue  problem 

[K]  -  [M]  =  [0]  (4) 

to  determine  the  mass-loaded  natural  frequencies  and  mode 
shapes . 

Purpose 

The  purpose  of  this  investigation  is  to  obtain  data, 
and  compare  the  modal  data  obtained,  using  several  modal 
prediction  techniques.  The  techniques  chosen  for  this 
comparison  are  those  discussed  previously,  namely;  (1) 
Glenesk's  use  of  Whaley’s  algorithm  (Ref  3);  (2)  the 
finite  element  method}  and  (3)  the  method  suggested  by 
Briggs  and  Whaley  (Ref  1).  Since  a  basis  for  comparison 
is  needed,  the  results  obtained  from  modal  prediction 
software  developed  by  Brown  (Ref  2)  will  be  used  as  a 
datum  in  error  precentage  calculations  for  methods  2 
and  3  presented  later  in  this  report.  However,  it  is 
felt  by  the  author  that  to  recalculate  error  values  for 
Method  1  which  would  be  based  on  a  different  datum  would 
be  unfair  to  both  Glenesk  and  Whaley's  algorithm.  Thus, 
all  values  presented  in  reference  to  the  use  of  Whaley’s 
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algorithm  will  be  directly  extracted  fjoin  Reference  3. 
The  same  complex,  rib-stiffened  panel  and  several  of 
the  discrete  mass-loading  configurations  investigated  by 
Glenesk  (Ref  3)  will  be  used  in  this  comparison. 

Ob  iectives 

The  objectives  of  this  investigation  are: 

(1)  Experimentally  measure  the  natural  fre- 
guencies  and  mode  shapes  of  a  complex  test  panel  in 
the  freguency  range  from  0-500  hz. 

(2)  Construct  finite  element  models  of  this 
panel  and  conduct  a  modal  analysis  on  both  unloaded  and 
mass-loaded  configurations. 

(3)  Use  equations  1,  2,  and  4  along  with 
unloaded  experimental  data  to  obtain  mass-loaded  nat¬ 
ural  frequencies  and  mode  shapes  for  the  test  panel. 

(4)  Present  a  comparison  of  the  results  of 


Objectives  1  through  3 . 


II  Modal  Analysis  and  To. si  iM'oi  •  's  i  'os 

Tost  Item 

The  test  item  was  a  panel  fabricated  usinq  draw¬ 
ings  of  an  upper  fuselage  panel  of  a  C-140  aircraft 
(Figure  1).  The  curved  panel  consisted  of  the  follow¬ 
ing  components : 

1 .  An  outer  skin 

2.  Five  longerons  of  two  different  cross  sections 

3.  Two  curved  m.ain  frame  ring  segments 

4.  Four  edge  doublers,  and 

5.  Various  attachment  hardware  and  bonding  to 
maintain  structural  integrity. 

Mass-Loading  Configurations 

Based  on  the  worst  case  errors  presented  by  Glenesk 
(Ref  3)  for  the  mass-loaded  panel,  four  test  configur¬ 
ations  were  chosen  for  comparison  of  the  three  methods. 
These  were  the  unloaded  panel  and  Glenesk *s  mass- 
loaded  configurations  2 ,  6,  and  7  (Table  1  and  Figure  2). 
The  unloaded  panel  was  included  as  a  mieans  of  compar¬ 
ing  the  change  in  mode  shape  with  natural  frequency 
which  occurred  between  the  unloaded  case  and  each 
mass-loaded  configuration.  Plots  of  unloaded  versus 
mass-loaded  mode  shapes  allowed  visualization  of  this 
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Panel,  USAF  Drawing  X704933. 


Table  1 


J  oad )  !u‘i 


Configuration 

>!ass  (lb) 

Xq  (ft) 

■SHB| 

2 

0.2920 

1.168 

1  .  168 

6 

0.9158 

1.667 

1.799 

7 

0.1199 

0.999 

0.915 

Fiaure  2.  Mass-loadinq  locations  (after  Glcnesk 


r-' 


0.279 


BAY  1 


■  -•  BAY  A 


Notes ; 

1. 

Dimensions  in  feet. 

2. 

/A  specifies  . -onf iquration  7 
tion. 

mass- loading  loca- 

3. 

BAY  1  specifies  internal  un 

.suppnrted  bay. 

chanqo  in  mode  sliape  (Figuin,'  i). 

Structural  Models 

Three  grid  sets  were  chosen  to  model  this  panel . 

The  first  grid  set  was  utilized  in  an  effort  to  verify 
Glenesk's  results  by  using  Glenesk’s  grid  set  and 
modal  analysis  procedure.  Glenesk's  model  did  not 
consider  the  discrete  components  as  separate  members. 
Instead,  it  accounted  for  the  total  panel  mass  and 
smeared  this  mass  over  the  grid  which  was  inset  some¬ 
what  from  the  panel  edges.  The  result  was  a  homogene¬ 
ous,  constant  thickness  flat  plate  with  25  grid  points. 
This  model  will  be  referred  to  as  the  Smeared  Stiffener 
Uniform  Model  (Figure  4 ) . 

The  second  grid  was  chosen  to  coincide  with  a 
finite  element  model  which  accurately  modelled  the 
discrete  structural  components  by  allowing  for  panel 
curvature,  discrete  member  cross-sectional  geometry, 
and  offsets  of  component  neutral  axes.  This  model  also 
consisted  of  25  grid  points.  It  will  bo  referred  to 
as  the  Discrete  Stiffener  Model  (Figure  5). 

The  final  model  was  chosen  to  investigate  the 
modes  of  vibration  of  each  of  the  internal  bays  (Figure 
2  ) .  This  grid  was  necessary  because  the  previous 
two  models  basically  ignored  the  motion  of  the  unsupported 
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internal  bays.  This  model,  which  con;  i  :ed  of  65 
grid  points,  will  be  referred  to  as  the  Bay  Modes 
Model  (Figure  6).  This  model  contained  the  previ¬ 
ous  two  models  as  subsets . 


Modal  Analysis  Test  Procedures 


The  Smeared  Stiffener  Uniform  Model  was  tested 


first.  The  Modal  Assurance  Criterion  (>L4C)  developed 
by  Brown  (Ref  3)  along  with  discrete  Transfer  Function 
data  were  used  as  the  basis  for  identifying  candidate 
frequency  ranges  from  which  the  natural  frequencies 
for  each  configuration  were  determined.  The  MAC 
function  is  defined  as  (Ref  3) 


MAC 


Syr-(()j)  ^ 

S  (  (jL*  ^  S  yy  (  (t)  } 


where  Syj-  =  the  stable  average  of  the  cross  power 
spectrum  between  two  response  measure¬ 
ment  points 

Sj-j.  =  the  stable  average  of  the  auto  power 
spectrum  of  the  stationary  accelero¬ 
meter  response 

Syy  =  the  stable  average  of  the  auto  power 
spectrum  of  the  moveable  accelero¬ 
meter  response. 

Note  that  the  MAC  Function  differs  from  the  more  commonly 
used  Coherence  Function  in  that  the  two  measurements 
in  question  in  the  former  are  two  responses  to  an  im¬ 
pulse  excitation  whereas  like  quantities  for  the  latter 


mmr 


BPi 


r 


would  be  an  impulse  excitation  input  loP  the  forced  re¬ 
sponse.  In  both  cases,  the  existence  of  a  mode  is  indi¬ 
cated  by  a  region  of  closely  spaced  frequencies  where  the 
MAC  or  Coherence  Function  is  essentially  equal  to  one. 

Since  the  MAC  and  Transfer  Function  data  were  essen¬ 
tially  identical  with  Glenesk's  results,  these  data  were 
not  reduced  to  obtain  natural  frequencies  and  mode  shapes. 
Instead,  Glenesk's  results  (Ref  3)  will  be  used  in  the 
techniques  comparison.  Table  2  contains  a  summary  of  the 
pertinent  data.  Sample  MAC  and  Transfer  Function  plots 
and  data  are  contained  in  Figures  7,  8,  9,  and  10. 

The  Bay  Modes  Model  was  tested  next  using  modal 
analysis  software  developed  by  Brown  (Ref  2).  This  soft¬ 
ware  allowed  the  user  to  select  any  subset  of  the  model 
being  tested  and  consider  only  the  data  relative  to  that 
subset.  Thus,  it  was  not  necessary  to  repeat  this  test 
for  either  the  Sm.eared  Stiffener  Uniform  Model  or  the 
Discrete  Stiffener  Model.  Representative  mode  shape 
data  obtained  for  these  models  are  contained  in  Figure 
11. 

Modal  Analysis  Data  Reduction 

To  determine  the  natural  frequencies  of  a  given 
model  and  configuration,  one  grid  point  was  selected 
which  was  believed  to  contain  all  of  the  modes  in  the 
frequency  range  from  0-500  hz.  That  is,  it  was  believed 
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Table  2.  Experimental  Natural  Frequencies  Obtained  Using  Modal  Assurance 
Criterion  (MAC)  Function  (Extracted  from  Ref  3). 


Figure  7.  Sample  Modal  Assurance  Criteria  (MAC)  Plot. 
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that.  th('  point  did  not  lii‘  on  or  r.c  .n  ■  .i  nodo  line  Coy 
any  mode.  Examination  of  the  real  and  imaqinary  parts 
of  the  Tran.sfor  Eunotion  for  that  arid  point  revealed 
the  natural  frt'quenries  for  that  particular  model  cind 
conf iqurati on .  Detailed  modal  analysis  aiid  test  pro¬ 
cedures  are  contained  in  .Appendix  B.  Data  acquisi¬ 
tion  proqrams  for  use  in  the?  HowloLt-Packrj rd  HP-a  ll  IB 
Fourier  Analyser  arc  presented  in  Appendix  C. 

Once  the  natural  frequencies  for  a  particular  n;odei/ 
conf  igura  t  ion  had  been  identir'ied,  the  Transfer  Function 
data  from  each  point  on  the  grid  in  quo.stion  were  reduced 
to  yield  mode  shape  vectors  for  each  mode.  For  purposes 
of  comparison  with  the  finite  element  method  the  discrete 
mode  shape  vectors  were  converted  to  a  formiat  identi¬ 
cal  to  NASTR.AM  ou‘‘put.  X'ASTR.AN  data  were  processed  using 
a  standard  graphics  package,  GCSXAST  (Ref  6)  ,  which  wa.s 
used  to  display  plots  of  the  undeformed  versus  deformed 
mode  shapes  in  both  unloaded  and  m.ass-loadod  conf  iqur>it  ioi 

Resul ts 

Modal  Assurance  Criteria  Function  tost  results  are 
presented  in  Table  2.  Cor rospondina  results  from  the 
Modal  Analysis  testing  are  shown  in  Tablc>  j.  Coniparison 
of  the  data  in  these  two  tables  ri'Vfmls  numerous  are, is 
of  disagreement  in  the  natural  frequencies  of  each  con¬ 
figuration.  This  apparent  discrepancy  can  be  explained 
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when  one  considers  the  tost  tor-hn!qLn> 
set  of  data  were  obtained.  The  MAC  tt.'sting  used  response 
data  from  a  fixed  reference  acceleroineter  and  a  moveable 
accelerometer,  and  15  impulse  excitations  which  were 
randomly  spaced  over  the  entire  panel.  The  spatial 
randomness  of  the  excitations  enhanced  the  probability 
that  every  mode  in  the  structure  would  be  excited  in 
that  not  every  excitation  would  lie  on  a  node  line.  In 
contrast,  the  Modal  Analysis  testing  used  a  fixed  ex¬ 
citation  point  in  conjunction  with  a  moveable  accelero¬ 
meter  to  measure  the  structural  forced  response.  Inherent 
in  the  latter  technique  is  the  assum.pt ion  that  the  chosen 
excitation  point  never  lies  on  or  near  to  a  node  lino 
and  thus  the  mode  shapes  obtained  from,  this  techniqus 
represent  a  unique  set  of  modal  data.  Therefore,  if 
the  chosen  excitation  point  lies  on  or  near  to  a  node 
line,  some  modes  may  be  "missed"  during  examination  of 
Transfer  Function  data  because  these  modes  were  never 
excited  to  begin  with.  In  this  respect  it  would  appear 
that  the  MAC  function  data  may  be  the  more  accurate  data. 
Since  a  comparison  of  modal  prediction  techniques  is  the 
thrust  of  this  report  this  comparison  will  be  carried  out 
only  on  those  data  for  which  corresponding  results  in 
the  MAC  Function  data  are  available.  This  comparison 
method  will  be  used  throughout  the  remainder  of  this 
report. 
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Ill  Tho  Whalov  Algor  i  I !  ip.  * 


Overview 

The  algorithm  developed  by  Whaley  (Ref  13)  and  sub¬ 
sequently  used  by  Glenesk.  (Ref  3)  substitutes  suitable 
expressions  for  the  kinetic  and  potential  energy  of  a 
flat  plate  into  Lagrange's  equations  of  motion.  Then, 
taking  the  first  variation  of  the  expression  for  the 
virtual  work  of  the  applied  inertial  loads  with  respect 
to  the  generalized  coordinate,  the  expression  for  the 
generalized  force,  Qj^,  is  obtained.  Tho  final  form  of 
this  expression  is 


6P 
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)  = —  =-M  d.  (x  ,  y  )  Z  0.(x  ,  y  ) — r 
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-M  R^“(x^,  y  )  Z  i— ^(x  ,  y  ) - , 
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d-^g 


(5) 


dt‘ 


Equation  5  contains  the  effects  of  the  inertial  forces 


due  to  the  added  lumped  mass  (Mq?  Rx>  Rv>  Xq,  and  Yq) 

which  are  proportional  to  the  second  derivative  of  the 

generalized  coordinate,  - ^  or  q.;.  These  effects,  in 

dt'^  -* 

turn,  become  additions  to  tho  mass  matrix  when  the 


eigenvalue  problem  is  solved  for  natural  frequencies 


and  mode  shapes . 
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Resul t  s 

The  results  from  Whaley's  algorithm  are  summarized 
from  Glenesk.  (Ref  3)  in  Tables  4,  3,  and  6.  The  data 
generated  using  the  MAC  Function  wore  merged  with  Whaley's 
algorithm  (Ref  13)  to  yield  the  Whaley  Algorithm  results. 
Glenesk  (Ref  3)  noted  that  a  comparison  of  unloaded  versus 
mass-loaded  mode  shapes  was  the  moans  of  determining  the 
unloaded/predicted  mass-loaded  freguency  pairings. 

Examination  of  the  data  in  Tables  4,  5,  and  6  reveals 
the  largest  percentage  error  between  actual  and  predicted 
values  occurred  on  Configuration  7  (Table  6,  Mode  1)  where¬ 
as  the  smallest  percentage  error  value  occurred  on  Con¬ 
figuration  2  (Table  4,  Mode  9).  With  the  exception  of 
Mode  1,  Configuration  7,  all  predicted  frequency  values 
were  within  +  7.5%  of  the  experimentally  measured  values. 

Table  7  presents  a  comparison  of  the  absolute  per¬ 
centage  in  frequency  shift  from  the  unloaded  panel  to 
each  of  the  three  mass-loaded  configurations.  The  largest 
percentage  frequency  shifts  occur  when  the  largest  mass 
was  located  on  an  unsupported  portion  of  the  panel  skin. 
With  the  exception  of  modes  8  and  9  for  Configuration  6, 
all  frequency  shifts  were  within  +  3'V,  of  the  unloaded  fre¬ 
quency. 
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Table  4.  Data  Results  -  Configuration  i  extracted  from  f 


Testina  Techniauo 


Percent 


Mode 


1 


Modal  Assurance  Criteria 


( 

174.08 

187.62 

204.00 

225.31 

243.32 
263.80 
282.45 
363.08 


Table  5.  Data  Results  -  Configuration  6  (Extracted  from  Ref 


Mode 


1 


Testina  Techniauos 

Modal  Assurance  Criteria 

Khalov  Alqorithm 

169.71 

166.22 

180.72 

175.34 

199.30 

191.23 

208.27 

205.63 

238.22 

235.20 

252.20 

251.18 

279.95 

278.44 

292.20 

271.14 

360.20 

338.70 

Percen 

error 


-2.06 

-2.98 

-4.05 

-1.27 

-1.27 

-0.40 

-0.54 

-7.21 

-6.11 
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Table  6.  Data  Results  -  Configuration  ,  (Extracted  from  R<‘f 

3). 


Mode 

Testing  Tochniguos 

Whalov  Alaorithm 

143.10 

167.77 

17.24 

2 

175.09 

180.20 

2.91 

3 

184.83 

191.16 

3.44 

4 

204.30 

206.98 

1.31 

5 

238.22 

237.33 

-0.37 

6 

257.50 

250.45 

-2.73 

279.88 

278.70 

-0.42 

292.76 

292.48 

-0.09 

361.71 

364.39 

0.74 
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Table  7.  Percentage  of  Frequency  Shift  Using  Whaley's  Algorithm  versus  Unloaded 
Measured  Data  (Extracted  from  Ref  3). 
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IV  Tho  Finite  Element  '  i 


Overview 

Several  excellent  finite  element  codes  for  struc¬ 
tural  analysis  are  presently  in  existence.  One  such 
widely  used  program,  NASTR^MJ  (A’Asa  STRuctural  Analysis, 
Ref  9,11),  was  selected  for  use  in  the  finite  element 
modelling  and  modal  analysis  of  the  structure.  This 
code  includes  the  general  20  degree-of-f reedom  guadri- 
lateral  elements  (CQUAD2)  and  12  degree-of-f reedom  bar 
elements  (CBAR)  of  which  the  test  structure  was  con¬ 
structed.  NASTRAN  also  contains  provisions  to  allow  for 
the  offset  of  the  neutral  axes  of  che  bar  elements  from 
the  grid  points  which  were  defined  at  the  midsurf ace  of 
the  panel  skin.  Thus,  the  cross-sectional  and  spatial 
properties  of  each  component  stiffener  could  be  included 
in  the  analysis.  Onlv  the  out-of-plane  component  of  the 
vibration  (i.e.  -  radial  component  for  curved  models 
and  z-component  for  the  flat  model)  was  investigated. 

Finite  Element  Models 

Three  finite  models  were  constructed  to  coincide 
with  the  test  grids  described  in  Section  II.  The  Smear¬ 
ed  Stiffener  Uniform  Model  considered  the  panel  to  bo 
a  flat  plate  (Ref  3).  This  model  did  not  extend  to  the 
panel  extremities  but  was  inset  somewhat  from  the  panel 
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edges  (Figure  2)  to  coincide  v/ith  the  grid  used  by 
Glenesk  (Ref  3) .  The  mass  of  the  entire  panel  includ¬ 
ing  stiffeners  was  smeared  over  this  grid  to  provide 
a  uniform  thickness  model  with  homogeneous  material 
properties.  Although  the  mass  of  this  model  was  identi¬ 
cal  to  the  overall  structure  mass,  this  model  was  more 
dense  than  had  it  been  extended  to  the  geometric  test 
panel  boundaries.  It  consisted  of  25  grid  points  from 
which  data  were  obtained  and  16  quadrilateral  elements. 

The  Discrete  Stiffener  Model  (Figure  4)  was  con¬ 
structed  to  allow  for  panel  curvature,  discrete  stiff¬ 
ener  geometries  (i.e.  -  different  cross-sections),  and 
the  offset  of  the  stiffener  neutral  axes  from  the  panel 
surface.  This  model  was  designed  to  faithfully  repre¬ 
sent  the  panel  from  a  structural  standpoint  while  main¬ 
taining  the  25  grid  points  of  the  Smeared  Stiffener 
Uniform  Model.  Forty- four  bar  elements,  25  grid  points, 
and  16  quadrilateral  elements  were  used  in  this  model. 

The  Bay  Modes  Model  (Figure  5)  was  included  to  in¬ 
vestigate  the  motion  of  each  of  the  four  internal  bays. 
This  model  was  an  extension  of  the  Discrete  Stiffener 
Model  in  that  two  extra  sets  of  five  grid  points  per 
bay  were  added  to  that  model  to  derive  t’nis  model.  This 
model  was  constructed  of  65  grid  points,  48  quadrilateral 
elements,  and  76  bar  elements. 

Appendix  C  contains  a  sample  of  the  Bulk  Data  Docks 
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which  woro  usod  to  qenorati'  srod  il  d  c  :dr  lhi"-  i‘  inode]  , 

Results 

Only  the  Bay  Modes  Model  re'sult.'-,  ('I'ablf^:-  8  throuqh 
11)  will  be  discussed  since  this  model  yielded  the  be  t 
representation  of  the  overall  panel  motion.  As  in  the 
use  of  the  Whaley  Algorithm,  a  pairing  of  unloaded  \-('rsu.s 
mass-loaded  mode  shapes  (Figure  12)  ’.vas  used  to  track  the 
change  in  natural  frequency  between  the  unloaded  panel 
and  each  mass-loaded  configuration.  Upon  observation  of 
Figure  12  it  is  seen  that  this  process  is  somewhat  subjec¬ 
tive  in  nature  in  that  the  mode  shapes  do  not  remain  com¬ 
pletely  unaltered.  It  is  left  to  the  discretion  of  the 
engineer  to  properly  select  the  mode  pairings,  and  herein 
lies  a  potential  source  of  error.  In  an  attempt  to  have 
this  mode  pairing  as  unbiased  as  possible,  the  author  con¬ 
sulted  another  engineer  to  independently  aid  him  in  this 
process.  Mode  pairings  which  were  not  in  agreement  between 
the  author  and  the  other  engineer  wore  discussed  and  a 
consensus  of  opinion  arrived  at.  The  results  from  this 
mode  shape  pairing  exercise  are  presented  in  Table  12. 

Here  absolute  changes  in  natural  frequency  from  the  un¬ 
loaded  to  mass-loaded  panel  range  from  essentially  .eero-',, 

(Mode  20,  Configuration  2)  to  -20. 3-".  (Mode  16,  Configura¬ 
tion  6).  Out  of  57  such  pairings,  the  frequency  shift 
from  the  unloaded  to  the  mass-loadi'd  panel  in  was  within 

I 

i  7.5-%  of  the  unloaded  panel  in  70  percent  of  tlie  pairinns  i 
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Table  8.  N’ASTKAN’  Modal  Data  -  Unloaded  Pant? 


Table  9.  .‘JASTRAN  Modal  Data  -  Conf  iqural  i  o!i 


Table  10.  NASTRAN  Modal  Data  -  Conf iqurat ion  6 


9B 


Figure  12.  Mod<'  ,  Bay  Modes  •■  iniVo  I- 1  (uvoiit  M 

for  A.)  eio.m  Fanol  (20G.  'G  Hr, )  ,  B.)  Ho 
ration  2  C'OH.ll  Hr),  h.)  oonf  iaural  uni 
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X'  The  rso  o  ~  rs('udM  i  ri'.'t ' r:  .■  J_n  1  hi '  _  )\'<'r\'  of 

the  Discrtho  >!ass  .  S  t:  l  :  f  ness  ,  and  ..<.a.iDinq 
Matrices  and  SoJ  ut  i  on  of  f  i'.o 
Mass- Loaded  Problem 


Overviov 

The  method  suggested  by  Briggs  and  Whaley  (Ref  1) 
solves  the  general  matrix-vector  differential  equation 
using  generalized  coordinates  to  obtain  the  definitions 
of  equations  1  through  3.  Then,  using  modal  data  ob¬ 
tained  from  standard  vibration  tost  methods,  and  suit¬ 
able  matrix  manipulations  (See  Appendix  A) ,  equations 
1  through  3  are  solved  for  the  mass,  stiffness,  and 
damping  matrices.  The  usual  method  of  solution  re¬ 
quires  a  square  modal  matrix  which  is  inv'orted  in  the 
solution  for  these  matrices.  However,  in  general, 
this  matrix  may  be  rectangular  with  more  rows  (grid 
points)  than  columns  (mode  vectors).  Thus,  one  would 
like  to  }iave  a  means  to  solve  equations  1  through  3 
when  the  modal  matrix  is  non-square.  The  pseudoinverse 
(Ref  8)  presents  such  a  method  for  inverting  rectangu¬ 
lar  matrices  and  thus  potentially  for  a  more  general 
solution  to  these  equations. 

Data  Reduction 

A  computer  program  was  developed  to  lake  tlie  ex¬ 
perimentally  determined  unloaded  modal  data,  compute 
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the  unloaded  modal  matrix  and  requir  -d  r.nverses,  and 
solve  for  the  unloaded  mass,  stiffness,  and  damping 
matrices.  Since  only  the  natural  frequencies  and  their 
resulting  mode  shapes  were  of  interest  in  this  study, 
the  standard  eigenvalue  problem  posed  by  equation  4  was 
solved  with  suitable  additions  of  mass  in  the  mass  matrix 
for  the  mass-loaded  natural  frequencies  and  mode  shapes. 
Appendix  B  contains  a  more  detailed  explanation  of  the 
data  reduction  process;  Appendix  C  contains  a  copy  of 
the  program  used  for  data  reduction. 

Results 

An  extremely  simplified  experimental  set  up  consist¬ 
ing  of  nine  grid  points  was  used  for  each  mass-loaded  con¬ 
figuration  (Figure  13)  to  demonstrate  the  validity  of  the 
computer  program  before  extension  to  the  m.ore  general 
pseudoinverse  case  ivas  attempted.  The  structure  was  test¬ 
ed  in  both  unloaded  and  three  discrete  mass-loaded  con¬ 
figurations  from  which  nine  f requency/mode  shape  pairs 
were  identified  for  each  unloaded/mass-loaded  configura¬ 
tion  (Figure  14).  The  nine  modal  vectors  were  used  to 
form  a  square  9X9  modal  matrix;  the  nine  natural  fre¬ 
quencies  wort?  used  to  form  the  matrix  c;i  the  right  side 
of  equation  2.  Then,  using  the  technique  described  in 
Appendix  B,  the  unloaded  mass  and  stiffness  matrices  were 
recovered,  and  a  quantity  of  mass  equal  to  the  added  mass 
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Figure  14.  Nine  po i n i  ronrse  grid  niodo  stinpc  ;,  i~7TnTTaTTT 
tion  7,  ''odo  1,  r nr  A.)  I'nloadi'd  p.inc]  flC)). 
Hz),  H,  )  lo.ided  P.inel  (13'','>77  Hz),  e.)  :qn,. 

Modal  '■^airix  FTedirlion  (]r)M,,M.17  Hz)  .itui  H.  ) 
Pseudfi  i  nv- ■  PTf'diot.  inn  (Lt)'‘,.n7  H;-,).  (I)!',' 

cd  patK'l  denf)(  ed  by  dot  Ic'd  linr'^,) 
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was  added  to  tho  ri,5  term  (i.o.,  t}’.<-  iss  loading  poir.i. 
See  Figure  13)  of  the  mass  matrix.  Finally,  equation  4 
was  solved  usinq  a  generalized  eigenvalue  routine,  ETCZ" 
(Ref  3),  for  the  mass-loaded  natural  frequencies  and  mod' 
shapes.  These  predicted  mode  shapes  were  then  compared 
to  the  actual  measured  mode  shapes  to  verify  that  tho 
program  had  successfully  predicted  both  the  correct  natu¬ 
ral  frequency  and  its  corresponding  mode  shape.  From 
this  mode  shape  comparison  (Tables  13,  14,  and  15)  it 
was  determined  that  the  program  using  a  square  miodal  ma¬ 
trix  had  correctly  predicted  the  frequency  within  ±  5~°o 
error  24  out  of  27  times,  and  with  ±  10%  error  26  out 
of  27  times.  These  data  are  presented  in  Tables  13,  IT, 
and  15. 

Next,  the  last  column  of  the  modal  matrix  and  the 
last  row  and  column  of  the  matrix  of  natural  frequencies 
squared  were  deleted  to  simulate  a  case  where  fewer  than 
"n"  modes  were  measured  where  "n"  is  the  number  of  grid 
points  (or  rows  in  tho  modal  matrix).  The  modified  9  X 
8  modal  matrix  was  then  inverted,  and  9X9  mass  and 
stiffness  matrices  were  qenc^rated.  Tho  mass  matrix  Wris 
perturbed  by  addilion  of  quantity  of  mass  equal  to  tb.e 
mass  loading  con f  i ■  ura t ion  at  the  5,5  location,  and  tlii^ 
general  oigenvalui?  p;'oblem  was  again  solved  using  EiGZl  . 
Upon  oxaminat  i.fjn  of  tlv'  actual  unloaded  and  predicted 
mass-loadt'd  pseudei  ir.e:-so  frequency  data. 


It  was 


found 


that  those  data  are  identical  to  six  d  cimal  places. 

Of  even  more  concern  is  the  fact  that  the  predicted 
mode  shapes  v>?ere  quite  different  from  the  measured  mode 
shapes.  Thus,  from  this  preliminary  investigation,  it 
did  not  appear  that  the  method  employing  pseudoinverses 
would  yield  valid  results  and  further  attempts  at  its  use 
were  abandoned. 

However,  it  was  demonstrated  that  the  method  using 
square  modal  matrices  will  yield  valid  results.  One  draw¬ 
back  in  the  use  of  square  modal  matrices  in  this  method 
is  that  when  one  increases  the  dimension  of  the  modal 
matrix  (i.e.  -  the  number  of  grid  points  or  rows),  a 
corresponding  number  of  mode  shape  vectors  and  natural 
frequencies  must  be  generated.  At  the  outset  of  this  re¬ 
port  it  was  poin<.od  out  that  for  larger  numbers  of  grid 
points  this  may  be  both  undesirable  and  that  it  may  not 
be  possible  to  obtain  a  large  number  of  natural  frequen¬ 
cies.  Thus,  there  is  a  need  to  limit  the  number  of  grid 
points  w'hen  one  models  the  structure  in  question.  Un¬ 
fortunately,  there  is  at  present  no  general  method  to 
predict  the  number  of  natural  f roquencics/modc's  within 
a  given  frequency  range.  The  solutif)n  to  this  problem 
may  be  to  initially  bt'gin  wildi  <i  simplified  nine'  point 
grid,  test  the  unloaded  structure  .uid  reduc('  the  Tran:;- 
fer  Function  data  to  obtain  thf'  number  of  natural  f re'- 
quencies/modes  in  the  specific  frc'guency  range  of  interest  . 
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This  number  of  natural  frequoncLos/r;  will  give  an 

indication  of  tho  maximum  number  of  grid  points  avail¬ 
able  if  one  wishes  to  refine  this  grid  and  use  the  square 
modal  matrices.  Since  tho  Transfer  Function  data  thoor(?t- 
ically  contain  the  same  information  for  any  grid  point 
unless  the  response  accelerometer  for  a  particular  grid 
point  was  situated  on  a  node  line,  it  is  not  required 
to  initially  increase  the  number  of  grid  points  above 
nine.  During  data  reduction  to  determine  tho  maximum 
number  of  natural  frequencies  (or  grid  points),  it  will 
thus  bo  advantageous  to  reduce  the  data  from  more  than 
one  grid  point  to  insure  that  no  modes  were  "missed"  duo 
to  a  given  accelerometer  being  inadvertently  placed  on 
a  discrete  node  line. 


dR 


VT  Discussion  of  the  >‘t'i  iif 


General 

As  with  any  endeavour  of  this  type,  a  learning  curv<> 
is  associated  with  one’s  ability  to  use  unfamiliar  equip¬ 
ment,  software,  and  testing/modelling  techniques.  This 
was  especially  true  for  the  author  who  had  no  previous 
experience  in  modal  analysis  testing,  with  the  associated 
data  reduction  techniques,  nor  \>fith  finite  element  modi^l- 
ling  procedures.  Thus,  the  author  was  in  a  unique  posi¬ 
tion  to  bo  able  to  evaluate  the  throe  modal  analysis  ti''c:i::i 
quos  which  are  the  subject  of  this  report.  Admittedly 
this  assessment  of  the  three  methods  will  be  only  one 
person’s  viewpoint  and,  as  such,  is  somewhat  subjective! 
in  nature.  Howo/er,  it  represents  the  viewooint  of  orie 
who  was  previously  uninitiated  in  this  area  of  expertise. 
The  evaluation  of  each  method  will  bo  presented  separate¬ 
ly  in  the  successive  paragraphs  of  this  section. 

At  the  outset  of  this  investigation  it  was  cited 
that  the  Modal  Analysis  Software  di^v^elopt^d  by  Brown  (Ref 
2)  would  be  used  as  a  datum  for  comparing  thf‘  method;-. 
However,  because  the  mode  sliapes  f 'r^m  Brown  ’  s  software 
and  NASTRAM  wore  quite  different,  such  a  comparison 
between  experimental  data  and  MASTRAM  data  was  not  pos;;i- 
blc.  Furthermore,  since  the  method  u.^;  i  ng  pseudo  i  in’i  r  :;  ' 
was  abandonod  in  favor  of  square,  in\-ertib]('  0  X  R  modal 
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matrices,  and  (.'xtension  to  citlicr  J.  or  6a-  arid  point.-: 
was  not  possible  due  to  a  lack  of  a  sufficient  number  of 
natural  frequencies  in  the  specific  range  from  0-500  Hz, 
a  comparison  of  this  method  to  Browin' s  softw-are  Vv'as  like¬ 
wise  not  possible.  Neither  was  it  possible  to  tie  Glonosk's 
modal  data  to  corresponding  modal  data  generated  by  Brown's 
software  due  to  lack  of  sufficient  modal  data  in  Glenesk’s 
report  (Ref  2).  Thus,  each  method  will  be  discussed  as 
fully  as  possible  while  keeping  those  limitations  in  ir.ind. 

The  Whalev  Algorithm 

The  results  obtained  from  use  of  the  Whaley  Algoritlim 
were  extracted  directly  from  Reference  3  and  are  presented 
in  Tables  4  through  6.  Because  the  author  did  not  direct¬ 
ly  use  this  algorithm,  the  author  is  unable  to  present 
an  evaluation  of  the  difficulties  encountered  in  obtain¬ 
ing  results  from  this  method.  How'ever,  when  one  considers 
the  modes  which  were  predicted  by  the  algorithm  v/hon  modal 
testing  did  not  reveal  such  a  mode,  one  can  envision  seri¬ 
ous  limitations  to  its  use.  As  previously  noled,  Gleivsk 
reported  a  degradation  in  algoritlim  pm-fornumce  as  tlic 
size  of  the  added  lumped  mass  incroasad  ndativaa  to  t!.s 
overall  structure  r.i'-:-:.  Thu:-,  an  as  yet  undet('>Tm' ns'd 
upper  bound  exisi:-;,  beyond  whirdi  the  predicted  mode  :-,ii>ipex 
and  natural  f requi  nf'ies  will  ceas«'  to  have  an  accept  .Hole 
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degree  of  accuracy. 


one  to  obi  ilr. 


Tho  USD  of  iliis  a2.goritln:.  requ 
modal  data  on  tho  unloaded  structure  before  the  mass- 
loaded  modal  quantities  can  be  predicted.  This  rc'quiros 
modal  testing  facilities,  tost  hardv.-aro ,  and  tho  nccessa 
software  for  conducting  tho  required  Fourier  analysis 
on  tho  Transfer  Function  dat<;i.  This  can  mean  a  signi¬ 
ficant  commitment  in  terms  of  resources  (manpower,  money 
testing  facility,  etc).  How'over,  -with  the  advent  of 
portable  modal  analysis  equipment,  th.is  investment  in 
resources  is  expected  to  decline  rapidly.  One  factor 
not  as  easily  evaluated  is  the  level  of  expertise  re¬ 
quired  to  successfully  conduct  the  required  modal  test¬ 
ing  and  subsequent  data  reduction  in  order  to  obtain 
valid  results.  It  was  tho  author's  experience  that  one 
must  rely  heavily  upon  "experts"  in  this  field  to  avoid 
the  pitfalls  of  questionable  tost  procedures  or  of  an 
invalid  data  analysis.  The  knowledge  required  to  suc¬ 
cessfully  conduct  tho  testing  and  data  analysis  is  one 
of  those  intangible  factors  upon  which  it  is  extromLc^ly 
diffii.'ult  to  place  an  evaluation. 

Thf?  one  seri(5us  limitation  inhert-nt  in  tho  usi  ■  of 
this  alqoril  hr!  is  t  hat  it  will  not  pr'^dict  m.ass-loaded 
mode  ;-;hape  ;.  ItiO  'ad,  it  considers  th.!‘  Ciodf^  shapes  tc^ 
be  utuiltf'red  by  I  !■■'  addition  cjf  mass  ie  the  structuri'. 
Although  unaltiT.  d  r!:od('  shrapes  wt're  ebierved,  tlu'y  u('re 


ratidoti'.ly  sc.i  1 1  on’d  thronnhcjui  L-iu  .  ^  rirnc'iii.^il  dat-t. 

It  is  the  opinion  of  the  author  upon  review  of  the  i:ux- 
perimental  data  that  oven  with  a  froquc'ncy  shift  the 
mode  shapes  were  altered  far  more  offon  tlian  not.  This 
is  not  believed  to  be  a  function  of  experimental  ir.otliodf)- 
logy  or  questionable  data,  but  rather  an  actual  occurr.'ncs^ 
in  nature. 

The  Finite  Element  >!othod 

Again,  as  in  the  case  of  the  modal  analysis  testing, 
a  learning  curve  was  associated  v;itli  the  use  of  the 
finite  element  method  as  presented  in  NASTR.W .  Unlike 
modal  analysis  testing  v.'here  tlie  methodology  used  to 
obtain  and  reduce  the  data  is  all  important,  the  author 
found  the  structural  modelling  technique,  especially 
modelling  of  the  boundary  conditions,  to  be  critical  in 
obtaining  good  results.  Since  the  structure  v;as  support¬ 
ed  by  bungy  cords,  the  support  conditions  fell  somr'whoro 
between  the  "free-free"  case  and  the  "clamped-clar;pod" 
case  at  the  panel  boundaries.  Thus,  it  was  not  possible 
to  clamp  the  finite  element  model  at  the  edoes ,  noi'  was 
it  possible  to  allow  the  pan(?l  to  b(?  totally  uniass  t. r  :  innd . 
Several  methods  of  m'-^delling  the  boundary  conditions 
were  used  to  detin'minis  the  mod'.'l  whiadi  most  clnsr'ly  ap¬ 
proximated  the  actu  il  boundary  c..; sJit  ions.  Among  Llu' 

were  t!w?  restraint  (Tf  rigid  btady  mt'dir^  by 


methods  tried 


niodall-i.nq  tlv  ■  bunqy  cords  Icji'q  ■  :  :f  r  rods  v.'ith 

stiffness,  the  use  of  tho  SUJ’ORT  card  or  SPCl  cards  to 
restrain  sui  table  deqree.s-uf- f rc'odoi::  at  throe  non-col  ine  .r 
grid  points  on  tho  panel  edges,  and  the  use  of  tho  SUPOltT 
card  to  restrain  all  six  deqroc'-of- f rf'odom  at  ono  select¬ 
ed  grid  point.  In  tlio  latter  case  a  COMMd  card  cas  used 
with  small  m.ass  and  inertia  values  at  tho  soloctod  grid 
point  to  alleviate  the  singularity  in  the  mass  matrix  im¬ 
posed  by  restraining  six  degreos-of-f reedom  at  the  same 
grid  poirit.  The  teciinique  using  SPCl  cards  vas  chosen  as 
the  best  method  to  model  the  supports. 

In  general,  each  finite  element  model  generated  a 
greater  number  of  natural  frequencies  (Tables  8  through 
12)  than  the  MAC  function  test  results  (Tcibles  2,  -1,  5, 
and  6)  indicated  wore  present  in  tho  frequency  range  from 
0-500  Hz.  How'ever,  when  the  finite  oloniont  results  (Table 
12)  wore  compared  to  Modal  Aiialysis  test  results  (Table 
7),  approximately  tho  same  number  of  natural  frequencies 
were  obtained  using  these  two  methods.  Although  the  quani- 
ty  of  frequencies  obtained  using  the  finite  elem.ent  method 
was  approximately  the  same  as  those  obtained  from  Modal 
Analysis  results  and  tho  frequencies  wore  numerically 
similar,  it  becrime  apparr^nt  upon  comparison  of  tlie  mode 
shapes  from  numerically  similar  frequencies  that  more 
numerical  similarity  of  the  natural  f ri'quencies  was  in¬ 
adequate  by  itself  as  a  criterion  for  cf^niparing  th^ 
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Modal  Analysis  tost  results  with  si!:.il  ■ results  f ron 
NASTRAN .  As  the  desired  moans  of  comparison  of  these  too 
sets  of  data  vas  a  "pairing"  of  "like"  mode  shapes,  i.t 
was  unfoitunate  that  tliis  was  not  possible  due  to  the 
large  discrepancy  in  mode  shape  versus  natural,  frequency 
when  the  two  sets  of  data  were  compared.  Thus,  a  bettor 
means  of  comparing  analytical  to  experimental  data  is 
needed  and  should  be  the  subject  of  further  research. 

This  failure  in  ability  to  compare  expGrim.ental  to 
analytical  data  presents  a  perplexing  problem  for  tl'ic  de¬ 
signer  who  wishes  to  minimize  the  vibrational  effects  on 
the  hardware  to  be  installed  by  placing  the  associated 
attaciiment  hardware  at  points  of  minimum  vibration.  Which 
mode  shapes  does  the  designer  believe?  Even  more  distres¬ 
sing  are  the  "additional"  modes  which  wore  computed  by 
NASTRAN.  Do  these  modes  really  exist?  Are  they  heavily 
damped  modes  which  exist  not  far  from  the  noise  floor? 

Were  these  modes  missed  in  the  experimental  testing  due 
to  an  inadequate  sampling  bandwidth?  Or  are  these  "addition¬ 
al"  modes  pur'-'ly  "s;^nithotic"  modes  generated  by  the  solution 
software  in  NASTR.AN  but  not  actually  existing  in  nature? 

These  questions  remain  unanswered  and  are  seen  as  a  severe 
limitation  of  tlio  "state-of-the-art"  in  vibration  analysis. 

Althougli  a  pairing  of  experimental  versus  analytical 
(finite  element,  NASTR.AN-generated)  mode  shapes  was  not  suc¬ 
cessful,  no  .such  lack  of  success  was  encountered  when 
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pairing  either  XASTR  AN- gone  rated  un]  e  'i  data  to  NASTR.'.',- 
generatod  mass-loaded  data  (Table  12  and  Figure  12)  or 
when  corresponding  Modal  Analysis  Software  data  (Table  i) 
were  compared.  Thus,  both  methods  appear  to  bo  internally 
consistent  when  data  generated  by  a  particular  method  are 
compared . 

Since  the  thrust  of  this  investigation  has  been  the 
comparison  of  techniguos  for  predicting  mass-loaded  natu¬ 
ral  freguencies  and  mode  shapes,  and  not  to  argue  the  vali¬ 
dity  of  the  experimental  baseline  data  versus  NA3TRAN  data, 
one  would  be  inclined  to  agree  theit  NASTR.AN  presents  the 
easiest  technique  for  the  experienced  designer  to  predict 
mass-loaded  modal  data.  The  basis  for  this  is  that  if 
all  the  NASTRAN-predicted  modal  quantities  really  do  exist, 
the  designer  has  at  his  disposal  all  such  data  without  the 
need  for  testing  which  might  load  to  incomplete  results. 

One  area  which  has  not  boon  explored  is  the  fact  that 
the  Modal  Analysis  Software  uses  a  least-squares  algorithm 
to  fit  the  Transfer  Functions  for  residut's  and  phase  angles 
in  a  given  frequency  range.  In  the  generation  of  the  Modal 
Analysis  Software  mode  shap<_v:,  any  phase  angles  •jcncrati'd 
due  to  a  frequency  shift  b(''tw('>en  Tixiruifer  Functions  were 
ignored.  This,  by  itself,  will  result  in  non-orthoqona I 
modal  vectors.  Since  NASTRAN  qen<'rates  a  set  of  mutually 
orthogonal  modal  vectors,  it  may  ce  possible  to  con'.pare 


the  Modal  Annlv<-;is  rial  a  to  tlie  NASTRAN  data  if  the  erf(V't 


of  phase  angle  (i.e.  -  frequency  shLfi  .vLth  grid  point 
within  a  given  mode)  is  included  in  the  display  of  the 
Modal  Analysis  mode  shapes.  This  capability  is  not  pre¬ 
sently  available  on  the  HP5451B  Fourier  Analyser  used  in 
the  experimental  phase  of  tliis  effort.  Should  this  capa¬ 
bility  be  added  in  the  future,  the  appart'nt  discrepancy 
between  the  Modal  Analysis  and  MASTR.AM  results  may  be 
eliminated . 

The  percentage  in  absolute  froqacncy  shift  from  the 
unloaded  panel  to  the  mass-loaded  data  are  presented  in 
Table  12.  Unlike  Whaley’s  Algorithjn  where'  all  but  two 
cases  were  within  +  3%  of  the  unloaded  panel  results, 

17  out  of  58  (29.3a.)  were  greater  than  this  +  3?i  frequency 
shift.  Even  more  disturbing  is  the  wide  range  of  froquen" 
cy  shift  indicated  by  mass-loading  for  some  conf igriration;: 
(refer  to  Table  12,  Modes  12,  16,  and  22).  Thus,  it  would 
appear  that  at  least  part  of  the  NASTR.W  data  may  be  ques¬ 
tionable  . 

A  considerable  expenditure  of  the  autlior's  time  and 
computer  resources  was  made  in  obtaining  these  results. 

For  comparison  purposes,  the  author  kept  records  of  the 
computer  resources  used  to  obtain  the  modal  data  for  each 
of  the  three  finite  element  models  as  v;ell  as  for  designing, 
debugging,  and  running  the  pseudoinvnrse  program.  These 
data  are  presented  in  Table  16.  It  should  bo  empliasized 
that  the  author  was  totally  inexperienced  in  the  use  of 
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NASTRAN  at  the  outset  of  this  of  fort;  rus,  the  cxpe't;d.  sure 
of  resources  for  the  finite  element  data  using  NASTR.W 
can  be  expected  to  be  reduced  by  up  to  75"u  by  an  engineer 
who  is  more  experienced  in  the  use  of  XASTRAX. 

In  view  of  the  fact  that  the  finite  element  method 
results  did  not  compare  favorably  with  the  experimental  re¬ 
sults,  one  must  answer  the  question  whether  the  time  and 
effort  invested  in  obtaining  these  data  was  commensurate 
with  comparable  investm.ents  from  the  other  two  methods 
vis-a-vis  their  respective  results.  The  answer  lies  in 
the  accuracy  required  by  the  structural  modification,  and 
by  the  vibration  sensitivity  of  the  hardware  to  be  install¬ 
ed.  In  the  case  of  electro-optical  devices  which  are 
known  to  be  inherently  vibration  sensitive,  it  v;ould  ap¬ 
pear  that,  time  permitting,  this  investment  of  resources 
would  be  justified.  However,  one  must  weigh  the  desirc'd/ 
required  results  with  the  method  used  in  each  individual 
case.  Thus,  one  of  the  other  two  methods  (Whaley  Algorithn 
or  Matrix  Recovery)  may,  with  other  less  vibration  sensi¬ 
tive  installations,  bo  optimal  when  all  factors  are  cvun- 


sidered . 


Tabic  16.  Computer  nosources  Used 


)tain  Uodt 


Parameter 

Finite  Elrr.r  nt 
Models^ 

Pseudo in verse 
Matrix 
Recovery^ 

Central  Pro¬ 
cessor  Unit 
Time  (sec) 

Total 

Average 
Per  Run 

Total 

Average 

Per  Run 

37590 

178.16 

3660 

19.26 

Input/Output 
Time  (sec) 

75155 

356.18 

4640 

24.56 

Computer  Cost 

$4461.46 

$21.14 

$280.92 

$  1.49 

These  data  include  post-processor  time  using  GCSXA3T  for 

models  up  to  158  grid  points. 

2  ..... 

The  data  are  for  primarily  nine  grid  point  models  and 

include  post-processor  time  using  DISSPLA. 


The  Use  of  Pseudoinverses  In  the  Recovery  of  the  Discrete 


Matrices,  and  Solution  o 


the  Eigenvalue  Problem 

The  recovery  of  the  mass  and  stiffness  matrices  and 
solution  of  the  eigenvalue  problem  was  considered  by 
Whaley  (Ref  16)  for  a  two  degroe-of-f roedom  system  and  th 
case  of  a  square  modal  matrix.  When  one  extends  this  |jr 
cedure  to  the  case  of  a  non-square  mod.il  matri.x,  trie  pro 
blem  becomes  more  .■;mplicated  becausr'  the  mass  matrix 
will  in  general  be  'uLly  populated  and  thus  th('  problem 
becomes  a  qenerali.sed  eigenvalue  problem  (Ref  t). 
simple  example  for  a  three  dei'iree-of- f  reedrim.  s>':;1('m  is 
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presontod  in  Appendix  A.  There  it  is  liewn  that,  vhen 
only  two  mode  shapes  and  corresponding  natural  frequencic's 
are  considered,  the  resulting  mass-loaded  natural  fre¬ 
quency  and  mode  shape  do  not  closely  approximate  tlic 
analytical  solution.  Thus,  it  does  not  presently  appear 
that  this  method  will  yield  accurate  results  ft^r  the  case 
of  a  non-square  modal  matrix.  Further  research  into  this 
area  is  needed  as  an  algorithm  (technique)  for  predicting 
mass-loaded  natural  frequencies  and  shapes  for  the  general 
(non-square  modal  matrix)  case  would  bo  extremely  useful. 

Summary 

The  time  invested  in  each  method  (inciudinq  the  use 
of  Brown’s  software.  Ref  2)  was  approximately  equal. 
Although  the  author  does  not  have  such  data  on  Whalrn-’s 
method,  it  is  assumed  since  Glencsk  completed  an  independ¬ 
ent  study  effort  leading  to  a  Master’s  degree  th<it  the 
effort  expended  in  the  use  of  Whaley’s  Algorithm  would  be 
comparable  to  the  other  methods.  Eacli  method  has  both 
drawbacks  and  good  points,  and  the  use  of  any  one  method 
over  another  will  have  to  be  neasonc'd  Vv'ith  a  groat  amount 
of  engineering  jud^pnont. 


VII  RGcommond'i  t  ions 

In  retrospect,  it  appears  that  the  test  item  selected 
was  much  too  ambitious  for  the  comparison  of  the  three 
modal  analysis  techniques.  One  suggested  test  item  would 
be  a  uniform  thickness  flat  plate.  Many  experiments  have 
been  performed  on  such  an  item  and  analysis  of  the  uni¬ 
form  flat  plate  is  included  in  many  elementary  vibrations 
textbooks  (Ref  5).  Once  the  methods  have  been  tied  to¬ 
gether  using  this  simplified  model,  modifications  to  the 
structure  could  be  added  one  at  a  time  and  the  process 
would  be  repeated.  Thus  one  would  gain  faith  in  the  three 
methods  and,  at  the  same  time,  gain  valuable  experience 
concerning  the  interactions  between  the  various  struc¬ 
tural  components. 

Along  the  same  line  of  the  gradual  increase  in  struc¬ 
tural  complexity  is  the  incremental  increase  in  the  added 
point  masses.  It  is  envisioned  that  this  small  incre¬ 
mental  change  in  added  mass  would  significantly  aid  in 
tracking  f requoncy  shifts  during  the  met?iod  comparison 
phase  by  allowing  bettor  visual  comparison  of  the  various 
unloaded  and  ma.ss-lo.idod  mode  shapes. 

The  fact  that  ph  ise  angles  different  from  zero  w('re 
present  in  the  experimental  data  qiv(?s  rise  to  the  qut's- 
tion  "Were  those  real  modf's?"  Speci  f  ica J  ]  y ,  this  qucrd.ion 
is  precipitated  bi'caurse  a  mode  sh.ipe  which  e.m  be  virarilly 
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observed  is  inherently  real-valued  wl '  a  no  phase  angle 
between  the  peaks  in  the  Transfer  Functions  at  a  given 
frequency.  If,  in  fact,  a  complex  mode  did  exist,  a  com¬ 
plex  display  would  be  necessary  to  properly  observe  the 
behaviour  of  the  structure.  Further  research  on  the  sub¬ 
ject  of  complex- valued  mode  shapes  is  beyond  the  scope 
of  this  report  but  should  be  pursued  in  the  future. 

One  area  using  the  square/rectangular  modal  matrix 
which  has  not  as  yet  been  investigated  is  the  addition  of 
a  structure  such  as  a  tripod  v/hich  connects  three  or  more 
grid  points.  Hence,  off-diagonal  torm.s  are  generated  in 
the  mass  and  stiffness  matrices  which  further  complicate 
the  issue  because  the  subjc.v't  of  finite  element  modelling 
of  the  tripod  structure  (or  any  other  added  structure) 
enters  the  picture.  It  is  the  author’s  opinion  that  this 
would  be  an  interesting  problem  for  one  who  is  interested 
in  furthering  the  research  in  this  area. 
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Mathematical  Considerations  and  example  Problem: 
In  the  Calculation  -md  ll-e'  of  tr.e  F'-a  ■udoinvf?r  a: 

Mathemat  ical  Cons  i  der.'itions 

In  the 

pa  pc 

r  by'  Penrose  (Ref  8)  it  is  shovai 

th.at 

for  any  matrix  A 

there  is  one  and  only  f:ine  matri 

satisfying 

the  : 

our  conditions: 

(1) 

AA'*‘A  =  A 

(A- la) 

(2) 

A'^AA^  =  A+ 

(A- lb) 

(3) 

(AA+)*  =  A-A-*- 

(A-lc) 

(4) 

(A+A)*  =  A+A 

(A-ld) 

The  matrix 

A+,  c 

ailed  the  pseudoinverse  of  A,  is 

the 

matrix  such  that  for  the  inconsistent  sot  of  equations 
A  X  =  b,  the  solution  x  =  A"*"  B  represents  the  optim<.il 
solution  to  the  least  squares  problem  A  x  =  b, 

Nov,  consider  the  problem  (posed  by  Noble,  Ref  7) 

A  X  =  b  where  A  =  BC,  and  A,  B,  C  are,  respectively, 
m  X  n,  M  X  k,  and  k  x  n,  and  all  three  matrices  are  of 
rank  k,  then  the  solution  of  A  x  =  b  which  miniminos 

a)  the  sum  of  tiio  squares  of  the  residuals  r  ‘  r, 
where  r  =  b  -  A  x,  and 

_  T  - 

b)  the  sum  ■>:"  i.b.c^  squ.iros  of  th.('  unkno'.;ns  x  x, 
is  given  by^  x  =  A^  c,  w’n.ere 

A"*"  =  (CC^)~^  b\  (A- 2) 

When  one  ..iccompl  iahr.' as  ,t  dfa/ompos  i  t  ion  r:);'  the  matrix  A 


G-1 


in 


form  A 


■•■Ao  form  A  -  LI'  whi'i-c'  A  is  dti  i:i  x  :  ::ati'ix,  I,  is 

.  .  . 

an  m  x  m  niatrix,  and  U  is  an  ni  x  n  iratrix,  B,  P  ,  C, 

rp  '■J'’ 

and  C  in  equation  A-2  can  be  replaced  ty  I.,  L  ,  U,  and 


r.T 


Equation  A-2  tt.en  be-f-omos  (Ref  10) 


A+  = 

(uu'^’)  ^  (l'^l)  ^ 

l'^'  (A. 

Example  Problem 

1  . 

Calculation  of 

'■<’udoi  nvo r^-n 

Armed  with 

equation 

A- 3,  le*-  us 

cons i dor  a n  ex- 

ample  problem. 

Let 

the  matrix  A  be 

ropresontod  by 

1.0 

0.5 

A  = 

2.0 

1.0 

3.0 

0.0 

Followinq  the  procedure  of  Strang  (Ref  10) ,  the  follow¬ 
ing  procedure  is  used  to  decompose  A  into  L  and  U  factors. 
Using  an  elementary  matrix,  E,  to  accomplish  multipli¬ 
cation/addition  of  rows  in  matrix  A,  and  a  permuiation 
matrix,  P,  to  interchange  rows  to  achieve  tion-tero  pivot 
elements  in  U  as  needed, 

PEA  =  U,  and  (A- 4) 

A  =  (PE)“^  U  =  LU.  (A- 5) 

The  decompos ' tion  of  the  given  matrix  will  now 


be  calcu  j-  ,  -  cd  . 


- 

- 

1 

P 

— 

- 

1.0 

0.  ■■ 

1 

0 

0 

1 

1.0 

0 . 5 

1.0 

0 . 5 

A 

= 

2.0 

1.0 

- ► 

_ 

1 

0 

2.0 

1.0 

3 

!:a= 

0 . 0 

0 .  n 

3.0 

0.0 

-3 

0 

1 

3.0 

0.0 

0.0  - 

1 . 

1 

L 

- 

- 

- 

- 

1 

- 

Not irg 

tl'a  t 

the  SI 

'•cond 

row 

of 

LA 

contains 

ci 

nr  n- 

-ra'ro 

pivot  olomctit,  a  row  oxch.inqr'  with  tile  thi. rd  row  of 


EA  is  accomplisiied  using  a  pormutati  )s  matrix  as 

P  (EA)  = 


1 

0 

0 

0 

0 

1 

0 

1 

0 

L 

1.0  0.0 

0.0  0.0 

0.0  -1.0 


1.0  0.5 
0.0  -1.5 
0.0  0.0 


=  u. 


To  get  L,  augment  the  m  x  m  PE  matrix  with  an  m  x  m 
identity  matrix  on  the  right  side  and  perform  ele¬ 
mentary  row  operations  until  an  identity  matrix  is  ob¬ 
tained  in  place  of  the  original  position  of  PE .  The 
matrix  now  occupying  the  previous  location  of  the 


original  identity  matrix  is  the  inverse,  (PE) 

[peJi]  = 


.^-l 


1 

0 

0 

1 

0 

0 

1 

0 

0  1 

- 

3 

0 

1 

0 

1 

0 

- ► 

0 

0 

1  1 

- 

2 

1 

0 

0 

0 

1 

0 

1 

0  1 

1 

0 

0 

1 

0 

0 

0 

1 

0 

2 

0 

1 

. 

0 

0 

1 

3 

1 

0 

1 

3 

2 


0 

0 


Thus ,  A  =  LU  becomes 

A  =  LU  = 


— 

^  n 

r  1 

10  0 

1.0  0.5] 

1.0 

0.5 

2  0  1 

0.0  -1.5:= 

2.0 

1.0 

3  10 

0.0  o.o! 

L  -j 

3.0 

0.0 

According  to  Reference  14  if  the  matrix  U  con¬ 
tains  rows  consisting  of  all  zero  elements,  one  can 
delete  those  rows  in  U  and  the  corresponding  colurnns 
of  L  to  obtain  new  factors  L  and  U  such  that 

A  =  LU  (A- 6) 

Thus , 

A  =  LU  = 


"1 

0" 

[1.0 

0,5 

■]  .0 

0.  rr 

2 

3 

0 

1_ 

lo.o 

-1.5_ 

- 

2.0 

3.0 

1.0 

0.0 
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o  o 


'rho  psoudoitix'orsi^  of  A  bocoini's 


=  u'^  (UU^)  ^  (l'^l)  ^  l"*". 


Continuing  to  obtain  A  , 


(A- 7) 


A 


+ 


0.5 

-1.5 


1.0 

0.5 


1.0  0 .  o’ 

1.25  -0.75 

“l-l  3' 

12  3 

= 

0.5  -1.5 

-0.75  2.25 

3  1 

0  0  1 

0  0  15 

18  36  -30 


Once  the  matrix  a"*”  has  boon  calculated,  one  must  veri¬ 
fy  that  equations  A-la  and  A.-lb  hold.  Since  A  and  a"*"  contain 
only  real  numbers,  conditions  A-lc  and  /I- Id  need  not  bo 
checked.  In  general,  all  four  of  conditions  A-la  through 
A- Id  must  be  checked.  Chocking, 


aa'^a  = 


'1.0 

0.5' 

1 

45 

'  0 

0 

15’ 

'  i.n 

0 . 5 ' 

‘1 .0 

0 

.  5 

2.0 

1.0 

18 

36 

-30 

2 . 0 

1.0 

2 . 0 

1 

.0 

.3.0 

0.0_ 

J 

_3.0 

0.0 

.3.0 

0 

.c_ 

-  A. 


Simularly , 


+  + 

A  AA 


45 


45 


■  0 

0 

15' 

'1.0 

0.5 

1 

1 — 1 

o 

o 

18 

36 

-30 

2.0 

1.0 

13  3()  “30 

3.0 

0 . 0. 

'  J  > 

L 

■  0 

0 

15' 

18 

36 

-30 

- 

A  . 

Thus,  the  matrix  A  represents  the  pianidoi  nvf'rse 


t  he 


origina].  matrix  A. 

T  + 

If  we  take  B  =  A  and  compute  B  ,  we  f±nd 


+  1  °  +  T 

B^  =  ^  0  36  =  (A  ). 

15  -30 


If  we  let  A  =  U,  B  =  U  ,  equation  3-a  becomes 


U  MU  =  BMA  =  I. 


(A-8) 


Premultiplication  of  equation  A-8  by  B  and  postmultipli¬ 


cation  by  A  yields 


+  +  +  +  +  + 

B  BMAA  =  B  lA  =  B  A  . 


The  solution  to  equation  A-9  is 


M  =  b'^Ia'*' 


(A-9) 


(A-10) 


provided  the  consistency  condition 


bb‘^ia'*'a  =  I 


(A-11) 


holds.  The  reader  can  verify  equation  A-11  is  valid. 
Thus,  the  mass  matrix  becomes 


ihA*  =  Bth  =  i 


■  0 

18' 

■  0 

0 

15' 

0 

15 

36 

-30 

X 

45 

18 

36 

-30 

^2  648  1296 

-540  -1080 


-  540 
-1080 
1125 


■  0.16  0.32  -0.267 

=  0.32  0.61  -0.533 

-0.267  -0.533  0.556 


(A-12) 


Before  one  can  be  sure  ('quation  A-12  rf'prosonts  the  least 
squares  solution  to  the  mass  matrix,  otu'  final  check  must 
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be  made; 


B  M  A  i  I .  f  A- 1  3  ) 

The  reader  can  verify  that  equation  A-13  holds.  A  simihir 
procedure  is  followed  to  calculate  the  unloaded  stiffness  and 
damping  matrices  using  equations  2  and  3,  respectively. 

Example  Problem  2.  Calculation  of  Mass-]^oaded  Mass  and 
Stiffness  Matrices 

Following  the  procedure  cf  Example  Problem  A-1,  consider 
the  three  degree-of-freedom  system  of  Figure  A-1.  Let  = 

M2  =  M3  =  1  and  =  K2  =  K3  =  100.  Neglecting  damping  and 
the  forcing  functions,  the  equations  of  motion  become; 


‘1 

0 

0“ 

/  •  « 

■  200 

-100 

0  ■ 

0 

1 

0 

-2 

+ 

-100 

200 

-100 

U7  =  p 

0 

0 

1 

1^3 

0 

-100 

100 

(xjl  p) 

These  equations  are  the  three  degree-of-freedom  eigenvalue- 
problem  for  this  system.  Upon  solving  equation  A-ld  for  its 
eigenval^aes  (natural  frequencies)  and  eigenvectors  (mode 
shapes),  the  results  in  column  two  of  Table  A-1  are  obtained. 


Table  A-1.  Analytical  Data  From  Throe  Dogreo-of-Freedom 
Spr ing-Mass-Damper-System 


Parameter 

Mi=M2=M3=1 

Ki=K2=K3=100 

M3=M3=1,  M2=2 

Ki=K2=K 3=100 

Unnormalized 
Modal  Matrix 

f,J2=12.47 

fi)3=18.0ig 

.j;i=3.813 
(jj^  =  ll  .  84  3 
?, 1^=15. 648 
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fl(t) 


-^1 

X 

■ 

-^'2 

ms 

.  .3 

'-1 

8 

■ 

^2 

hid 

'-3 

bid 

Equations  of  Motion; 


Mass  Matrix  Damping  Matrix 


’M]^  o'  0 

Kl 

(C1+C2  “C2  0 

l''<l) 

0  M2  0 

x^ 

+ 

-C2  (C2+C3  -C3 

-■<2 

0  0  M3 

^3 

0  -C3  C3 

(■^'31 

Stiffness  Matrix  Forcing 


(Kt+K2) 

-K2 

0 

I'^i 

-5^2 

-K3 

l'’^2 

fa  0 

0 

-K3 

^3  J 

(^3 

f3(t) 

Function 


Figure  A-1.  Three  Degree  of  Freedom  System 


Next,  let  a  mass  of  magnitude  one  be  added  to  mass  M2  sucl 
that  M2  is  now  equal  to  two.  This  will  simulate  a  mass- 
loaded  configuration.  Data  for  this  configuration  arc  listed 
in  column  three  of  Table  A-1.  Note  that  the  mass  additiot. 
lowered  the  natural  frequencies  and  modified  the  mode  shapes. 
Now,  consider  the  first  two  natural  frequencies  and  cor¬ 
responding  mode  shapes  (Table  A-1,  column  2)  to  be  the  mi‘as- 
ured  modal  data  in  the  frequency  range  of  interest.  Using 
the  procedure  of  Example  Problem  A-1  to  recover  tht^  mass;  and 
stiffness  matrices,  the  "generalized  invert  equations  of 
motion  bf’come: 


0 


7 


r 


.485  - . 105  - .  13r 

->l| 

:._12  l^.vlO 

1 1 

-.105  .353  -.237 

^2  + 

-34.212  47.009  -22.523 

-.131  -.237  .249 

••1 

^3) 

_  11 .719  -22.523  12 . 871_ 

^3] 

|o| 

M 

After  adding  a  unit  of  mass  to  the  M22  term,  and  resolving 
the  eigenvalue  problem,  the  following  results  were  obtain¬ 
ed: 

=  0.02 
C1J2  =  4.45 
=  12.46 

Comparison  of  these  results  with  column  two  of  Table  A-1 
reveals  that,  with  the  exception  of  the  values  are 

nearly  identical  with  the  unloaded  results  .  and  that  the 
magnitude  of  the  frequency  shift  was  nowhere  near  that 
which  was  encountered  in  column  three  of  Table  A-1. 
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APPENDIX  B 


APPE.\'D1X  n 


Detailed  Modal  Analysis  and  Test  Proceduros 

Modal  Assurance  Criteria  (MAC)  Function  Testing 

The  tost  panel  was  suspended  from  a  ceiling-mount¬ 
ed  conduit  in  the  Structural  Vibrations  Branch  (Bldg  21C) 
of  the  Air  Force  Wright  Aeronautical  Laboratories  (FBG/ 
AFWAL)  by  bungy  cords  which  were  approximately  four  feet 
in  length.  These  bungy  cords  were  attached  to  the  panel 
through  holes  in  the  skin  located  near  each  of  the  four 
edges.  The  test  grid  used  by  Glonosk  (Ref  3)  v.'as  marked 
on  the  upper  surface,  and  a  reference  accelerom.eter  was 
mounted  on  the  skin  (see  "x"  on  Fig  2)  using  double-backed 
tape.  The  reference  accelerometer  and  moveable  accel¬ 
erometer  responses  were  individually  amplified  and  filt¬ 
ered  through  a  pair  of  matched  filters  before  their  signals 
were  input  to  the  HP-SdiilB  Fourier  Analyser  (Ref  ■!). 

The  reference  accolorometer  was  used  to  provide  the  trig¬ 
ger  signal  for  data  sampling. 

With  the  panel  configured  as  described,  I  biO  move¬ 
able  accelerometer  was  mounted  at  grid  point  numbi:u'  orn  • . 

The  panel  was  then  excited  (tapped)  13  times  at  randtmil  y 
selected  locations.  Following  each  tap,  and  prior  to 
computations  using  th('s<'  d.ita  in  the  I'ourJer  Analyser, 
each  data  sam.ple  was  viewed  ^.ep.tral  cl  y  on  ref  la'sh- 1  raci^ 
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oscilloscope  to  insure  the  tappinq  i;  d  :  ot  over loadc'd 
the  internal  electronics,  or  that  multiple  excitations 
had  not  occurred  during  the  data  sample.  Assuming  .m 
acceptable  data  sample  had  been  obtained,  these  data 
were  stored  in  the  Fourier  Analyser.  Eac?i  time  a  good 
data  sample  was  processed,  the  response  data  ior  that 
sample  was  averaged  with  the  previous  sample(s)  such  tliat 
after  15  iterations  the  stable  average  for  each  channel 
was  computed.  Using  these  average  response  data  with 
the  reference  accelerometer  as  the  input  and  the  move¬ 
able  accelerometer  as  the  system  forced  response,  the 
Fourier  Analyser  calculated  both  the  MAC  and  Transfer 
Functions  for  that  grid  point  and  displayed  the  Transfer 
Function  on  the  oscilloscope.  At  this  point  the  user 
had  the  option  to  either  accept  the  data,  or  reject  the 
data  and  reaccomplish  this  process  for  tliat  grid  point. 

If  accepted,  hard  copies,  including  plots  and  dictitial 
iniormation,  were  made  of  Transfer  Futiction  and  MAC 
Function  data.  The  Transfer  Function  data  were  stored 
for  later  use. 

The  next  stop  was  to  succossivoly  mount  t!m'  mo'.o- 
able  accclorr)mf!tor  at  each  of  tiie  rottuiininq  arid  points 
and  repeat  this  procedure,  Onc('  th<?  on;  ire  grid  h  id 
boon  survi’yr'd,  the  MAC  Function  data  wrue  ex, ii. lined  'o 
d<atermine  c.ind  id.itf'  frequf’ncif's  r  ry;-  v.Tiicii  a  curw'  i  1 
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of  the  Trans'*'!'  Function  data  v.'oiild  be  :  ic.c .  Prior  to 
curve-fitting  the  Transfer  Function  data,  each  discrete 
Transfer  Function  was  oxaniinod  at  the  given  fri^quency 
range  to  determine  the  data  with  the  "cleanest"  spike. 

This  spike  was  fit  first  using  a  J.east-squaros  algorittiiri 
to  determine  the  natural  frequency  and  damping  ratio  to 
be  used  in  each  of  remaining  Transfer  Functions  to  deter¬ 
mine  the  mode  shape  at  that  frequency.  This  process  was 
repeated  for  each  candidate  frequency  rang*!  identified  by 
by  the  MAC  Function  data  in  the  range.'  of  0-500  Hz. 

Modal  Analysis  Test  Procedure 

Several  grid  sets  were  used  in  this  testing:  the 
25  grid  point  scheme  of  Glenesk  (Ref  3,  Figure  2),  a  25  grid 
point  scheme  (Figure  5)  to  coincide  with  the  addition  of 
structural  components,  and  a  65  grid  point  model  with  fine 
meshes  superimposed  on  each  of  the  four  internal  bays  (Figure 
6)  to  investigate  the  modes  of  vibi'ation  of  each  of  the.si:' 
bays.  Since  tlie  tost  procedure  is  identical  for  each  grid 
sot,  it  will  only  be  discussed  fc^r  the  25  grid  point  s-^iiome. 

This  testing  used  a  mo\'e._ible  accelerome!  ('r  to  -;ui'\'cy 
the  grid.  However,  unlikf'  th*'  FAC  Tunc)  ion  t'';-t!ng,  a 
second  .stationary  accfG.eromet('r  w  i.;  not;  usf'd .  Instead, 
one  grid  point  was  sr-lf'cted  (numb<'r  one  foi-  this  grid) 
at  which  t.h*'  structure  would  br'  exritf'd  by  an  inpacl 
hammer  configured  with  a  load  cc'll  t  c.i  m<'asur('  1  h'-'  amplitude 
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of  the  impact.  The  output  fi  or;,  th,'  it,:  .  ur-jr-d  a-, 

the  trigger  source  witli  the  movoablt'  accolcromotor  meas¬ 
uring  the  forced  system  response.  Again,  unlil-io  the  .‘'.AC 
function  testing,  the  structure  was  exc.i  ted  only  five  time 
at  each  grid  point.  The  number  five  Vvs^s  chosen  for  two 
reasons;  first,  it  was  noticed  during  the  MAC  function 
testing  that  more  than  five  data  samples  had  no  signifi¬ 
cant  impact  on  increasing  the  accuracy  of  the  average  rSD; 
and  second,  due  to  the  largo  number  of  grid  points  used  in 
the  larger  models,  including  more  than  five  iterations  pc'r 
grid  point  would  cause  excessively  long  testing  while  pro¬ 
viding  little,  if  any,  improvement  in  the  quality  of  the 
data . 

The  process  began  by  mounting  the  accc-loromet(^r  at 
grid  point  number  one  and  exciting  the  structure  using 
the  impact  hammer  five  times  at  the  selected  exfitatiori 
point.  As  in  the  MAC  function  tosLing,  the  data  wore 
examined  following  each  discrete  sampled  to  determine 
if  an  overload  of  the  internal  eloctroiiics  or  a  multi. pli.' 
excitation  had  occurrt'd.  A.'-^.surninq  acceptable  dat.i 
obtained,  the  dirscreto  force  input  was  n;ult.iplied  a 
force  window  (Figure  B-1)  to  artificLaJJy  rc:)rce  it  to 
die  out  rapidly  as  wouJ.d  be  the  case  for  an  impulse  i'x~ 
citation.  similar];.,  I  hi'  structural  foi'cr'd  r('^;pons(' 
wa.s  multiplied  by  ati  I'xpc.ineiit. rally  di'cayinq  windo'n-  (!  iciuii’ 

B-2)  to  .simulate  the  effr'et  of  daP’-ping  et;  t  h('  ; ;  i  ru( -t  u  r,  1 1 
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response . 


The  a\'erage  of  ea'S'h  s  trip!' ■  ■„  '  h  pre-coodiiei 
samplf's  was  computed  artor  oacli  of  the  five  discrete  d..it,_; 
samples  was  accepted.  As  before,  .a  Transfer  huiiction  vers 
computed  for  each  point.  Upon  c’omplrh. ion  of  rhe  fifth 
iteration,  the  average  Transfer  Functi^on  and  a  Cohorf'nce 
Futiction  were  computed  and  individually  displayed.  A 
decision  wliether  to  accept  the  data  was  based  on  a  Coher¬ 
ence  Function  with  many  of  its  values  at  or  near  one  in 
the  0-500  Hz  range.  This  process  was  repeated  until  the 
entire  grid  had  been  sui-veyed. 

Once  the  survey  was  complete,  a  grid  point  was  select¬ 
ed  which  was  believed  to  contain  all  the  modes  of  the 
structure.  Several  methods  for  examining  the  Transfer 
Function  data  at  that  grid  point  were  available  in  the 
software.  Among  * hese  methods  were  a  Kenncdy-Plancu  CireJf 
fit,  a  loast-squ£U'os  algor itlnr,,  and  an  option  which  merely 
computed  the  magnitude  of  the  Transfer  Function  at  'eloctod 
f reg  (.’ncies  Vvhi'ro  the  real  and  i.nuigitiury  parts  of  the  Tiaans'- 
fer  Funr‘  ion  wore  00  degrees  out  of  piuise  with  eacri  olru-r. 
Tnitiatlly  the  author  .selected  the  i.'iagii i  tud.'  option  a;-  it 
ignores  dariping.  However,  tiae  least-squares  alar)!' i  t  h:",  v.  is 
even t  u a  1 1  y  u s e d  to  a n h  1  y.s o  t  h<^  d u  t  .u  . 


Mg  t  hi^mg  I:  i  r-  1 1  -  ft  \po  r  i  r'^en  t  ;  1  i’rora^d! ;  g,  > 

This  procedur*^  uses  tin'  expe  r  i  ru  ■[', !  a  !  1  y  d"!aM-minsd 
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Analysis  Section  presented  earlier  in  t})is  Appendix.  The 
mathematical  portion  of  this  procedure  begins  with  the 
general  matrix-vector  differential  equation  used  in  modal 
analysis , 

[M]  X  +  [C]  X  +  [K]  X  =  f(t),  (B-l) 

where  [m]  is  the  mass  matrix,  [C]  is  the  damping  matrix, 
and  [k]  is  the  stiffness  matrix,  all  of  which  are  square 
and  symmetric  matrices.  A  generalized  coordinates  solu¬ 
tion  to  this  equation  results  in  "n"  uncoupled  equations. 
Assuming  the  solution  to  equation  1  can  be  written  in  the 
form 

x(t)  =  [U]  q(t)  ,  (B-2) 

where  [U]  is  a  square  matrix  consisting  of  mode  shape 
vectors  as  its  columns  (i.e.-the  modal  matrix),  then  it 
follows  that  q(t)  is  the  generalized  coordinate.  Following 
the  derivation  of  Meirovich  (Ref  5)  after  substitution  of 
equation  B-2  into  equation  B-l  and  premultiplication  of 
equation  B-l  by  [U]"^,  the  following  definitions  are  made: 


Eu]'^[mj  [U]  =  [I] 
[U]'^[K]  [U]  =1  wj, 

[U]'^[C]  [U] 


(B-3a) 


(B-3b) 


(B-3c) 
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with  these  definitions  equation  B-1  becomes 


qj-  +  2  ^  r^r^r  "**  'X'r^r  ( B"4 ) 

rp 

where  fr(t)  =  [U]  f(t)  is  the  generalized  force. 


Equation  B-4  can  be  solved  for  at  most  "n"  natural 
frequencies.  However,  there  are  often  less  than  "n" 
natural  frequencies  in  the  frequency  range  of  interest 
such  that  the  modal  matrix  contains  fewer  columns  than  it 
contains  rows.  It  is  proposed  to  premultiply  each  of  the 
equations  B-3  by  and  postmultiply  each  of  these 

equations  by  [u]”^  in  such  a  way  as  to  isolate  the  mass, 
damping,  and  stiffness  matrices  on  the  left  side  of  equa¬ 
tions  B-3  which  involve  these  respective  matrices.  In 
general,  [U]"^  and  [U]  will  be  non-square  (rectangular) 
matrices  necessitating  the  use  of  the  pseudoinverse  devel¬ 
oped  by  Penrose  (Ref  8)  to  find  their  inverses.  Using 
a  superscript  "+"  to  indicate  a  pseudoinverse  and  after 
the  previously  mentioned  pre-  and  post-multiplications, 
equations  B-3  become 


[M]  [U]  [U]'^ 
[K]  [U]  [U]"" 

[C]  [U]  [U]"^ 


[I]  [U]'^  (B-5a) 
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Equations  B-5  can  be  further  reduced  usiu;  standard  linear 
algebra  matrix  manipulations  (Noble,  Ref  7)  to  yield  the 
mass,  stiffness,  and  damping  matrices  respectively  as 

[M]  =  [I]  [U]'"  {B-6a) 

tc]  = 

Reference  8  (Penrose)  contains  the  basic  theory  of  the 
pseudoinverse . 

With  these  definitions  for  the  mass,  damping,  and  stiff¬ 
ness  matrices,  and  the  experimentally  determined  modal  data 
for  the  unloaded  panel  gathered  in  the  Modal  Analysis  Section 
presented  earlier  in  this  report,  one  can  calculate  the 
mass,  damping,  and  stiffness  matrices.  Appendix  A  contains 
a  simple  example  of  the  technique. 

Since  the  thrust  of  this  procedure  is  to  accurately 
predict  the  mass-loaded  natural  frequencies  and  mode  shapes 
of  a  complex  structure  when  only  the  unloaded  modal  data 
are  known,  the  procedure  must  be  continued  to  recover  mass- 
loaded  modal  data  (Wn*"^  n»  Un) .  For  example,  consider  the 
mass  matrix.  When  the  designer  has  made  his  decision  re¬ 
garding  the  placement  of  the  additional  mass  (hardware)  in 
the  structure  in  question,  he  need  only  add  an  appropriate 


n  n 


:) 


[u]'" 


(B-6c)  . 


(B-6b) 
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mass  in  the  mass  matrix  at  the  proper  grid  location  to 
obtain  the  modified  mass  matrix.  A  similar  procedure  is 
followed  to  obtain  the  modified  stiffness  and  damping  ma¬ 
trices.  Assuming  one  is  only  interested  in  mode  shapes 
and  natural  frequencies,  the  modified  mass  and  stiffness 
matrices  are  substituted  into  equation  4.  Equation  4  is 
then  solved  for  the  mass-loaded  natural  frequencies  and  mode 
shapes . 


Finite  Element  Modelling  Procedure 

Finite  Element  Models.  The  three  finite  element  models 
used  in  this  analysis  were  described  in  Section  II.  Because 
the  development  of  a  finite  element  model  is  described  in 
Reference  11,  the  main  question  became  how  to  add  a  quanti¬ 
ty  of  mass  to  the  mass  matrix  which  would  be  used  to  solve 
for  the  modified  natural  frequencies  and  mode  shapes. 

NASTRAN  provides  this  capability  in  the  form  of  a  C0NM2 
card  which  adds  a  finite-valued  point  mass,  including  its 
inertia  properties  and  neutral  axes  offset,  to  a  particular 
location  in  the  mass  matrix.  Another  card,  the  ASETl  card, 
was  used  to  discard  all  but  the  out-of-plane  translation 
from  the  analysis  set  in  the  solution  for  mode  shapes  and 
natural  frequencies.  Appendix  C  contains  listings  of  a  sam¬ 
ple  NASTRAN  deck  used  in  the  modal  analysis  of  each  con¬ 
figuration. 

Analytical  Procedures.  NASTRAN  (Ref  11)  uses  several 


rigid  formats  to  provide  flexibility  to  the  user  in  the 
analysis  at  hand.  Rigid  Format  3  was  chosen  for  the  modal 
analysis  of  the  test  panel.  This  format  neglects  damping 
and  solves  equation  B-1  with  the  forcing  function  and  damp¬ 
ing  matrices  set  equal  to  zero.  The  output  from  this  For¬ 
mat  is  a  table  of  the  "n”  natural  frequencies  and  tables  of 
the  "n"  eigenvectors  (mode  shapes).  These  data  were  then 
post-processed  using  GCSNAST  to  be  displayed  on  a  standard 
computer  terminal . 
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FB'iTlOO,  tni00*c*<l‘‘0000.  AT  HI  KSCN, 


ATTACH.MSL  .in-LI«RAR  Y  ,SN-ASf>. 
iTTACHfOISSRl* ifP-L  fBH 
t|Bfl*RY,T-SL,nTSSPLA, 

FTN.PL-lOOaO. 

LCO. 

ATTACH,PLT10  3B,ID-A780?B'5**;N-ASn, 

PLT1038. 

REMIN0«TAPF91. 

POUTf  ,TAPF<?9,TIO-AF  •  P  I  D-F  B  I  %  S  T«C^  A  .OC -PR. 

*•  END  OF  RFC0»0 

P«OCRA'»  P4TRIX(  INPUT, nUTPUT#PtFItf-0.TAP£S-|NPUT,TlPE6»nuTPUn 

C#«« «#•••«•« 

c* 

CP 

CP 

CP  this  is  a  test  version  of  N4TRIX  USING  THE  NINE  CRIO  POINT  HOOEL  ANO  THE 
CP  resulting  "rof  shapes  ANO  natural  FREQUENCKS  to  verify  the  validity  of  THr 
CP  PROCRAN  TO  CALCULATE  BASS-tOAOEO  NATURAL  FRCOUFNCfES  ANf)  **OOE  SHAPfS,  NINE 
C*  separate  paSS-LOAOEO  cases  are  CONSIOERFO  in  ACCOROANCf  with  THF  cases 
CP  TESTfO  BY  GLENFSK^  SEVERAL  CASES  FOR  EACH  CONFIGURATION  ARE  CONSIOERED: 

CP  THAT  IS»  several  "STARTING  POINTS"  ARE  CONSIOfRED  FOR  THE  RANGE  OF  NATUR- 
CP  AL  FRFOUFNCIfS  TO  SfULATF  A  RFAL-MORLO  CASE  WHERE  ONE  IS  ONLY  INTERESTED 
CP  IN  A  SPECIFIC  FREQUENCY  RANGE. 

CP 

CP 

CP 

C44444444*4«44444«444««4444444«4444444444444«4444444444*44444444444*44444*4444 


CP 

CPTHIS  PROGRAB  SOLVES  THE  GENERAL  «ATR Ix-VECTOR  OIFFFRENTfAL  EQUATION 
CPFOR  THE  BASS,  OAWPINC,  AND  STIFFNESS  "ATRICES  USING  CENfRALUEO 
CPCOOROINATES.  THE  FOLLOWING  OEFlNllinNS  ARE  "ADEt 

CP 

CP  UT  •  AB  4  U  •  At 

CP  UT  P  AC  P  U  •  7ETA 

CP  UT  P  Aft  P  U  •  OBECA 

CP 


c» 

WHERE  U 

THE  NOOAL  BAIRIX 

c* 

UT 

a 

THE  TRANSPOSE  OF  THE  BCDAL  "ATRIX 

c» 

AB 

■ 

THE  NASS  NATRIX 

c» 

At 

• 

the  identity  BATRIX 

c« 

AC 

• 

THE  OABPfNC  BirfffX 

c* 

AR 

• 

the  stiffness  HATRIX 

c* 

zeta 

• 

A  diagonal  “ATRIX  with  THF  PRODUCT  OF  THE 

c* 

NATURAL  FRfOUfNCY  ANO  THE  OABRInC  RATIO  FOR 

c* 

EACH  BODE  OROEREO  IN  INCRfaSING  "ODE  NUMBER 

c* 

ALONG  THE  BAIN  DIAGONAL  AND  FEROES  ELSFwmERE 

c» 

OBECA 

• 

A  niACONAL  B4TRIX  WITH  THf  SQUARE  DF  THt 

c» 

NATURAL  FRCOUCNCY  FOR  EACH  BOOf  ALONG  THF 

c« 

"AIN  DIAGONAL  ORDERED  IN  INCRFASINC  BODE 

c« 

NU"BER  AND  ZEROES  ELSFWHfRf, 

CP 

C44444444444444444444444444444444444444444444444444444444444444444444 

CP 

CPSOLVING  THPSF  EQUATIONS  USING  STANDARD  "ATRlX  ••  A  N I  PUL  A  T  I  HNS  YIELDS: 


CP 

CP 

AB 

•  UCITN 

P  Al 

♦  UGIN 

CP 

AC 

•  UGITN 

P  ZETA  ♦  UGIN 

CP 

AK 

•  UCITN 

♦  n* 

'EGA  4  UGIN 

CP 

CP 

WHERF  UGIN 

•  THF  GFNFRALTZEO 

INVERSE  OF  Thf  NORBALIZEO 

CP 

NODAL  NATRIX 

CP 

UCITN 

•  THE  transpose  Of 

UGIN 

CP 

CPNOTE 

TO 
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IF 

J  IS  THF  NUNBER  of 

G«T0  points  ANO  L  IS 
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CP 
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CP 
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CP 
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CP 


Figure  C-1.  Program  Matrix  Listing, 
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c*N»o  •  Nij'»foirit  fiOM  oi«fN5inN 
•  NU*««F9  OF  BQmS 
C^NC  •  NU'**»F#  Cf  COLUMNS 
C*INlT!fcLTfF  THf  maTOUE^: 

Pl^FNSTON  uuf<J,«)i«oiFi^i*r)iFri<>>,oifii<»i*oir?i<)i 

01  SESSION  XAOt»YA«3l,/Af1l,X|<)i»Vf<9l«/Of9l«r)rC?>.rrf3l»Z/n» 

C1"€NSI0H  XX»|9|,YyYI<>l«;7?<‘»l 

Ot**FNSION  AALP<9| 

OI«*r^5  rON  U5f  9«9)  ,06  f9,<»»  ,07(9, 9>  ,0419,9) 

01 NS  ION  WN40AI 91 , WAtF19l «UA ( 9 , 9  I  ,( « R Oft  1 9 ) 

01  "f  NS  ION  AA  M  9,91  ,F  19,91,019,91  •  Hf  9 , 9  I  ,  0«<f  C  A  1 9 , 9  I  ,Uf9,9),ur(9,9) 
OTNFNSION  ACI9,9t,An9,9|,FETA<9,9l,FFf9,91,C0|9,g|  ,HH|9,9»  ,HA(9,9 
l),FYF|9,9|,CCC19,9|,Cf9,9l ,0 1 9 ,9 ) ,€ 1 9, 91 , / E T AS  1 9 ) .ONE  CAS 1 9  I , AL 1 9  I 
Pr^ENSrON  Eft*>0RU91 
CINENSION  WAMeoA19),llCITNt9,9t 
OMENS  ION  SU**19|  ,SU**l<9l,SU"?19t,SUN3<9l 
OMFnSION  SUNA19I 

REAL  ANN|9«9t  ,U3P1R,9|«SY14t«W(0116l,U319,4t»4f.  TA?f9t,WK4(l6?l 
ftEAL  Af9,9|,0ClN19,9l,TriL«Sf  9)  •Wftl)6^>«Anf  9,9),AK|9,9I  ,BbTAf9|  ,U?I 
19,9I,VX;IM2I  ,AAN19,9  t  •AA*(t9,9),RETA119|  ,WK1  I  l6^k,Ull9,9»,W"ECAf9) 
Cn-ftLE*  ALFA  191 , 7 1 9 • 9 1 « ALE  A t <9 1 , / 1 19 ,9  I ,ONf C AO  1 9 1 , Bt AN90 A ( 9 1 
CONPLEY  ALE A2f9|,/3|9,91,ALPI91 

OMENS  ION  ON  EC  It  8,81  ,n«<EGAlt  8t,Atl8,SI,CC19,8),U3TP|9,81,ANA19,9l, 
100(9,81 
NR-9 
NC»9 
ND«8 

CALL  CONftftS 
NCASES-3 

00  R  Kft-1,NCASES 

C9 

«EAO  structural  crio  points 

C9 

CAP* AAAA AAAAAAAA A* AAAA A* ••AAAAAAAAAPPAAAAAPAAP A  AAA* 

00  20?  Ml, NR 
ftEAOA,Y( I  I .Yl n , 70(1  I 
fE(F0F(5HNPUTI.NE,0l  STC® 

20?  PftfNTA,Y(l(,Y(M, 70(11 

CPAA«AAP«AAA«AA««AAAA««««««««««P«#«#P««A«PP««**PP*P**AAAA»AA*AA*A*AAAAAP 

CA 

€A  INITIALIZE  THE  NATRICES 

CA 

(♦PPPPAAAAAPPAAPAAAAAPAAAAAAAAAAPAPPAAAAAAAAAAPAAAAAAAAPAAAAAAAAAPAAAAAP 

00  111  Ml, nr 
00  111  J-1,NC 
All,  JM0,0 
U<  I,  JM0,0 
U2(l, Jf*0,0 
UA(1,  JM0«0 
UT(1,IM0«0 
UOINt  J,T  MO.O 
UGITNI  I,  iMO.O 
ZETASt  JMO.O 
ALf  JMl.O 

111  ONECASIJMO.O 
00  2  MlfNC 
DO  2  J-t,NC 
AKIflMC.O 
ONEGA!  I  •  JM0«0 
.?  2ETA(1,|MC,0 

00  333  Ml, NR 
00  333  J-l,NR 
ANf  I,  JM0«0 
AK1I,JM0.0 
333  ACdfJMO.O 
00  TO  Ml, NO 
00  TO  J-l,Nn 
AIIMJMO.O 
TO  0»  rCKI,  JMO.O 

PRTNTA,«TME  MITlALlZfO  MATRICES  ARE2** 

PRINTA,  •  THE  nOOAL  *«ATRTX:« 

DO  f»l  Ml, NR 

A1  Pt|NTA,fUfl«i»,J«l,NC1 


Figure  C-1,  Program  Matrix  Listing  (Continued) 


POINT***  TMF  *00*t  ooTOIlf  TOONSPOM:* 

no  h?  )-l.NC 

6?  PIIINT*,IUTI  J,ll.l-ltNO» 

POINT*,"  TMf  CfNFOALITfn  INVfOSt!" 

cn  n  1-1. Nc 

6J  PRINT* , tUCINI 1,1 1 , I -I ,N»» 

PRINT*,"  THf  CfNIOOLl/fO  IN»f  RSf  TRANSPORT:" 

00  NA  1-1. NR 

(,A  PR  1NT*  ,  lUCITNI  I  ,  II  .  1- 1  ,Nr  I 

PRINT*,"  TMt  ItTA  NATRlx:* 

00  I-l.NC 

PRINT*,I It TAII ,  II  I J-l .NCI 
PRINT*,"  THF  0"tGA  NATRIXt" 

no  hh  1-1, NC 

bt>  PR  INT*,I0NFCAI  I  ,  11  .  1-1  .NCI 

PRINT*,"  THf  •DU"NT*  "ATOICfS:" 

print*,"  ThF  AL  "ATRIX-." 
on  ni  i-i,NC 
111  PRINT*, Aim 

PRINT*,"  TMt  PA'i'i  NATRIXt" 

00  11?  I-1,NR 

11?  PR  I  NT*, I  API  I , II .  1-1 .NR  I 

POINT*,"  THF  n*"P|Nr,  NATOIX:" 

00  in  i-i,No 

113  po INT* , I  AC  I  I  .  1 1 .  J- 1 .NO  I 

PRINT*."  THF  STIFtNfi^  "ATRIX:" 

00  n*  1-1, NO 

11*  PO I  NT *, I AK I  I  . J1 ,  1- I  .NO  I 

PRINT*,*  THI  IOTNIITT  -ATOT*:" 

00  ITS  I-l.NC 

111  PR  INT*,  I  AM  I  .  II.  I- I  .NCI 


Fiquro  C-1.  rroqi-um  Matrix  Listiiiq  (Contimu'd)  . 
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PB|NT#.SU''f  J» 

JI-SORTISU«<  J»  I 
on  !-l,NR 

2\B  «(  ! »  Jl  -Af  I ,  Jl 

00  ^?l  I-l,NR 

?2l  PR INTA  •  f A| T  ,  Jl ,  1*1 ,NC I 
00  ?2n  I-l,N9 
00  220  J-UNC 
UI1«  M-AII,  Jl 
220  UT  t  J«n«At  I ,  il 
Caaaaaaaaaaaaaaaaaaaaaaa 
CP 

CP  OflfTF  THF  lA^T  CPLU*'N  FRO"  THf  lOOAt  "AfRiX  U  TO  f  OK  ••  THf  RECTANGUIAR 
C*  •MATRIX  U3. 

CP 

CAAAAAPAAAAAPAPAPAPPPPPPPPPPPPPPPPPPPP 

00  fll  !-ltNR 
00  «l  J-ltNO 
01  U1 I I«  Jl-At  I  •  Jl 

PR!NTPf"THlS  IS  THF  U  "ATRix:** 
no  Oil  I  -UNR 

911  PR  TNT* • (U( 1 « J I f J« 1 «NR I 

PRINT*, "THIS  IS  THf  A  MATRIX:" 
on  9\  I-1,NR 

91  PRINT*, ( Af  I  •  J  I  • J- 1  ,NR  I 
PflINTP,"THIS  IS  THf  UT  NATRIXJ" 
on  92  1*1 , NR 

92  PRINT*,! UT  U  «  Jl « i-l «NR I 
PRINT*, "THIS  IS  THf  U1  naTRIX:** 

00  02  I-1,NR 

02  PRlNT*,!Ul!l,JI,J-l,Nni 
00  71  I-l,NR 
71  ALP(l)-0,0 

on  03  i-i,No 
0"f Cl  1 1 , 1 l-ONf CAl ! I  1**2 
61  AlCI,II«t«0 

PR|NTP,"TH!S  IS  THf  0*»Fr.l  "ATRIX?" 

00  84  I-1,N0 

64  PRINT*, !rPFCl!NJI«J«l*Nni 

PRINT*, "THIS  IS  TMf  A1  NATRfX:" 

00  84  t-l«N0 

04  PRINT*, ! A1 n , Jl , J«1,N0I 

c******************************************************* *********** •*••***** 

c* 

c* 

C*C0«PUT6  THf  GCNfRALI/fO  INVfRSF  OF  THE  NOOAL  NATRiX  AND  ITS  TRANSPOSE 

c* 

c* 

c*************************************************************************** 

NRni-9 

NR0«9 

NC-9 

NR«9 

TOl-0.0 

20  CAlt  LCINf  I  A,NPO,NR,Nr  ,Tnt,UCIN,NRPl,S,)«K,If  Rl 
PRINT*, "THIS  IS  THf  -MATRIX  UClNZ" 

00  77  I-l,NC 

77  PR1NT*,(UC.tNI  I  ,  Jl  ,  J-l  ,NR  t 
00  1214  I-l,NB 
no  1214  J-l,NC 
1214  UCITNf I ,  JI-UCIN(  J,  I  I 

PRINT*, "THIS  IS  THf  xaTHJX  UCITN:" 
nn  414  I -1 ,NR 

414  PR  INT*,!IIGI  TNI  I  ,  Jl  ,  J-  I  ,NCI 

print*, "THf  VAlUfS  OF  ffTA  ANp  p-fGA  ARf,  Rf S P F C T | V f L Y : " 

00  14  i-i,Nr 

ONf CAS! I  I -P-f  GAS !  M / I  00* 

14  PR!NT*,/fTASII»,r*“^,AS!II 

c****************************** •****«***•*•*• **•«***•«* . . 

c* 

C*  CONRUTF  TMf  CFNFRAVIMP  INVFRSf  pf  THI  "ATRIX  in.  Tmjs  IS  THf  MATRIX  U1P, 

c* 

c***# ************************** ***•*♦*•••••****•••**•••*•******************• 

N*m«9 


Figure  C-1.  Proqr'ain  Matrix  Listing  (LotiVinued)  . 
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NCI-* 

TOL-O.O 

NRO*-8 

CALL  LC!NF !U3  *NH0  3,NR  3 .NCItTOL  *U3P,NttOA * S  3 * WK  3. I( » 1 ) 
PRINT*, "THIS  IS  IHF  -ftTRIU  U3P:" 

DO  TO  l-l,NO 

30  PR f NT*, IU3P ( I «  it •  i-1 tNC I 

€• 

C*  TRANSPnSf  f/3  TO  fORR  UlT.  ALSO  TPANSPOSI  U3P  TO  FORH  U3TP, 

CP 

C«*«*p*«««*«*«*«*«««*«**«««***«**«********#« *•»*••»•*•• **••*# 

00  31  1-1, NR 
00  31  J-UNO 

31  U3TP  I  I  ,  H-U3P(  J,  I  I 

PRINT*, "THIS  IS  Thf  "iTRIit  U3TP:* 

DO  3?  !-l,Nfi 

3?  PRINT*, IU3TPI 1 , Jl , J-l,NOI 


C*POR"  THP  7FTA,  ONFCA,  ANO  At  NATRtCfSi 

C* 

C* 

c******************************************* 

PRINT*, "THf  actual  "ATRlCfS  ARft" 

00  S  1-l.NC 

•PTAM,ri-?*/FrASII  t  •0"r/;AS  f  I  t 

PfCAt  t  ,  1  I-O-^ECAST  I  t*•^ 

Af(f,tl-1.0 
^  CONT'NUF 

PRINT*,"  THf  /fTA  NATRtf:" 

OC  2^  l-l»NC 

25  PR  INT*,(  Tf  TAIN  ( 1  I  ♦  M  •  I  «NCI 

PRINT*,"  THE  0*fCA  lATRfX:" 

00  1-1, NC 

26  PRINT*,  IO*»f  CAM  ,  1 1  t ,  I  I  •?  ,NCt 

PRINT*,"  THf  lOfNTITT  *ATR|i:" 

00  if  I-t,NC 

27  PRINT*,! At f I , it ,  i-l ,NC t 
00  Si  l-UNR 

00  At  IT-1,NC 
CU  .111-0,0 
00  Al  R-l.NC 

ccT,ni-cn,in*uciTN(i,*(t*AiTK,ni 
♦l  CONfTNUf 

00  A 2  I- I  .NR 
00  *2  ir-i,NR 
ANM,tlt-0.0 
00  *2  K-l,NC 

ANf  t  ,1  II  •1*11  ,  T  I  t^CI  I  .RI*UO|NIR,  Ml 
AANII  ,  n  )-A»»I  I  ,I  I  1 
A2  CONTINIJF 

no  A3  I-l ,NR 
00  A3  II-1,NC 

oiifiit-r.o 

00  A3  R-lfNC 

Of  I  ,111-011 ,1  I  l♦UCITN^  ,R  l♦7FTAf  K,H  • 
A3  CONTINUF 

DO  AA  t-l,NR 
00  AA  ll-lfNR 
AC(l,Ill-0.0 
00  AA  K-1,NC 

ATI  I  f  I  M  -  ACI  I  ,11  »  •on  ,R  l•l»G^N<R,  I II 

AA  CONTINUE 

00  A5  f-l,NR 
00  A5  I1-1*NC 
€f  I  ,f ri-0.0 
DO  A5  R-1,NC 

6  I  I  ,1  1 1  -  E I  I • n I *001 TNf ( ,* 1*0*6  GAIK, II  I 

A5  CONTINUE 

00  A6  1-1  ,NR 
on  A6  t1-l,NR 
ARC  I  tf 1 1-0,0 


Figure  C-l,  Program  Matrix  Listing  (Continued)  • 
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UC ITN( I |K (♦UCI Nl K  ,11 


00  *6  lOliXC 

46  CONTINUE 

00  I  I-l.NS 

on  I  j*i.NC 

IE(»«SI»N|I,JM.LE.10.E-8»  »4|1.J»-0.0 

lfl*BS(»KII,lll.LE.10.f-81  INI  1,11-0.0 

I  If (IBS  I »C I  I , n I .LC . lO.E-8*  »CII,J»-0.0 

PRINT*,”  the  D*«PINC  ”»tru;” 

00  48  I-l,HR 

48  PRINT*,  I  *C  H  ,  1 1  1  J”1  •Nl’ ' 

PRINT*,”  the  STIfFNESS  N*TRIx:” 

00  40  I-l.NR 

PMNT*’,”CHECn*C0HS1STEnCY  CONOITIQN  FOR  THE  IDENTITY  N4TBI*: 
00  000  I-l,NR 
00  OOO  J-1,NR 
PM,  JI-0.0 
DO  000  K-l,NC 

000  f  I  I, J I -El  1,1 1 
00  00?  I-1,NR 
00  00?  1-I,NC 
HI  I ,  1  1-0.0 
00  00?  K-l.NR 

00?  H(I,H-M(1,1I  ♦  f  I  I  ,K  I  *U(K  ,  1 1 

00  001  I-l.NC 
00  003  1-1, NC 
kkl II , n -0.0 
DO  001  N-1,NR 

003  *4  I  I  I  ,  1 1 - »* n  1 • 1 •  ♦  UT I  I ,NI*HI K , 1 1 

00  0003  l-l,NC 
00  0001  1-1, NC 

0003  If  I i RS 14 1 M I  ,  1 1  I . t E. 1 0. t-8 I  44111,11-0.0 
00  004  I-l,NC 

004  PRINT*, (44l(  I, II, 1-1, NCI 

PRINT*, ”CHfCN  CONSISTENCY  CONDITION  FOR  THE  ?fT*  H4TRI*- 
00  005  I-1,NC 
DO  005  l-l,NC 
FFII,  ll-O.O 
00  005  N-l,NR 

005  FF  I  I  ,  ll-FFI  I  ,11  ♦  UTI  I  .KI*UCITNII<,1I 
00  006  I-l,NC 
00  006  l•l,NC 
CCI I , II -0.0 
or  006  «-l,NC 

006  CCI  I  ,  II -CCI  I  .  j)  ♦  FFII  ,NI*?fT*lll,ll 
OP  OOT  I-l,NC 
OC  00?  I-1,NC 
HHI I, ll-O.O 
00  00?  K-l,NR 

007  HHI  I ,  I  I -HHI  I  , 1  I  »  UCINI  I  ,NI*UIK,  II 
00  008  I-1,NC 
00  008  I-l.NC 
H4 I  I ,  II -0.0 
00  008  N-1,NC 

008  H*l I ,  II -H*l I ,  II  ♦  CCI I ,KI*HHIN, II 
00  0004  I-l,NC 
00  0004  J-l,NC 


0904  |F(48SIH1II.1II.IE.]0.E-8I  H4II, JI-0.0 
00  000  I-I,NC 

000  PR INT*, IH4I I , 11 ,  I- l.NC I 

PRINT*. -CHECK  CONSISTENCY  CDNOITION  FOR  THE  OKEC*  "ITRIX:” 
DO  1003  I-I,NC 
00  1003  1-1, NC 
FFF II , 11-0.0 
00  1001  K-l.NC 

1003  FFFII. Jl-f f f I  I , I  I  ♦  f f I I,NI*nNEC*(K. 11 
DO  1001  1-1, NC 
00  1001  1-1, Nr 

CCCI I .11-0.0 
00  1001  K-l.NC 

1001  CCCI I ,  I  l-CCCI 1 , II  •  FF F I  I ,K l*HH|K .1 1 
00  0005  I-l.NC 
00  OOOS  1-I,NC 

0095  IF IIPSICCC I  I , II  I .L E . lO.E-K I  CCCII, JI-0.0 


Figure  C-1.  Program  Matrix  Listing  (Continued) 
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r 


—  'Tgfl'r!.\,Ng - - 

100?  PR|NT**rCCC( ( •  It *i«l «NCI 
PRINT*, "THf  NATRIv:* 

00  too*  {•UNR 
100*  PRINT*,(ANU • Jl, J«l,NCI 

C* •****#«**«*•**«*«*♦•«**•**•#•*♦•♦•♦*•***•*•*#*•••••*••*••••••**•••• •*•**•*•* 

c* 

c* 

c* 

€♦  SINCf  THf  NJSS  WA5  AOCEO  TO  CRfO  POINT  NUNBfk  FIVE  ON  THE  ACTUAL  PANEL,  THE 
C*  "ASS  NATRIW  NiLL  Pf  PERTURAfO  IN  THf  5,6  POSITION  Bt  AOniMC  A  QUANTITY  OF 
€♦  NASS  fOUAL  TO  THf  NASS  LOADING  ON  Thf  ACTUAL  PANEL.  THIS  kiLL  SINULATE  THE 
C*  "REAL  NORLD*  CASE.  NOTE  THAT  THE  NASS  WAS  PLACED  ON  ONLY  ONE  GRID  POINT.  IF 
C*  THE  NASS  NAS  A  STRUCTURE  SUCH  AS  A  TRIPOD  NHlCH  CO^NrCTEO  TO  THREE  GRID 
C*  POINTS,  THE  ADOEO  STRUCTURE  NOUL D  HAVE  TO  BE  "ODfLLEO  IN  THAT  IT  NOULO  AOD 
C*  BOTH  WASS  AND  STIFFNESS  TQ  THE  PANEL.  THIS  AOOE D  WASS  AND  STlffNfcSS  «OULO 
C*  GENERATE  AOCITIONS  TQ  OFF  DIAGONAL  TERNS  Of  LIKE  SU'^PfR  TO  THE  TwO  CONNfCT- 
C*  FO  GRID  POINTS.  AS  PREVIOUSLY  NOTED,  THIS  PROGRA*  SOLVES  THE  GENERALI/EO 
C*  EIGENVAIUE  PROBLEN.  THUS,  THERE  IS  NO  NEED  TO  CREATE  A  SIHULAR  "OOEL  FOR 
C*  ANT  ADDITIONAL  OA«P|NC  INCURREC  BY  THE  ADDITIONAL  STRUCTURE. 

C* 

c* 

c*********«***« •«••••«*•«•«•«**#♦*••*•«•«••••««•«•*•«•*••••••«•••*•#«•*#•#• A*# 

PR INT* ."ENTER  The  loop  TO  PFRTUWfl  THE  HASS  HATPI*:" 

(F( KR.EQ.lt  AH(S,6|<AAK(S,6l«6.1t4* 

IFIKK.EQ.?!  AHl6,6fAA^tS,St*0.?*20 
IffKK.fO.Tt  ANI^tSI^AAHIS.SMO.^ISB 
PRINT*, ■  THf  NASS  HATRIV2" 

00  SO  I»1,NR 

50  PRINT*,llAH(t, n ,J«1,NCI 

PRINT*, "THE  HOOIFIfD  HASS  HATRIK2" 

00  *7  I-l,NR 

47  PR INT*, < AHf 1 , it , j-l ,NR t 

PRINT*, -THE  STIFFNESS  "ATRIY:" 

00  5fl  I«1,NR 

5B  PR  INT*«  f  AKU  «  J  t ,  J«l,NC  1 

C**«*«********* *«•««« A  A************************* . A* A*** •♦••♦••♦♦•♦♦*** 

c* 

€•  SOLVE  THE  CEnERALWEO  EIGENVALUE  PROBLEM  FOR  THE  HASS-LOAOfO  NATURAL  FRE- 
CP  OUENCIES  and  "00£  shapes  NITH  The  square  NOOAL  natru  U. 

c* 

C***a**a««a*aaa***a*aaa**#aaaa«aaaa*aaaaaaaaa*aa>vaaa#aaaa*aaaa*aa****a*aa«aaaa 

1AN«R 

NANAK«<» 
t/«9 
I  JOB-? 

CALL  E  rC/FlAK, (AK.ANfl AN.NANAK ,I  JOB, AL F A , BE T A , 7 , 1  7 ,WK? ,I £R?> 

N-<J 

PRINT*, "THESE  ARE  THE  VALUES  OF  RLAHBOAS" 

00  l?34  l-l,N 

RLAH80AI  n-CSORT  IALFAI  1 1  FBETAI  1 1 1 
PRINT*, RLA4B0AI I t 
1?34  CONTINUE 

PRINT*, "THESE  ARE  THE  VALUES  OF  THE  7  NAi’RIX;" 

00  lOOD  t-l,N 
1009  PRINT*,<7<I, Jl, J-l,Nt 

PRINT*, "THE  PREOICTEO  "OOE  SHAPES  ANO  NATURAL  FREQUENCIES  FOR  AN|5 
l,5t-",ANT5,5t,"ARE:" 

00  ?l  l•l,NC 

?1  NAKBOAI I l•CSQRT^RLANBOAtII♦*?t 
PR  INT*  ,1  VAmBOAM  I  •  1«L,NCI 

LA-9 

CALL  VSRTAfWAwpOA.LAI 
PRlMT*,tWA«flOA(I f , |-l,NCf 
00  ??  t-l,NC 
00  ??  J*1,NC 

fFIABS  INAKPDAI 1 t-RFALI RLA-ROAf it tl.LE. lO.E-0 t  CO  TO  ?3 
IFIABSCWAKBOAU I-KEALIRLAN0OAI Jl I t.OT,10,E-B>  GO  TO  ?? 

23  00  ?4  K-l,NR 
U?(K, I t-C.O 

24  U?tK,M-U?(K,1l*REAL<7(K,JII 
22  CONTINUE 

00  15?  J-l,NC 
SUN4I JI-O.O 
00  153  !-l,NR 

15  3  SUN4I  J>«SU-4Ut  »U7M  ,  J  t  * •  7 


Figure  C-1.  Program  Matrix  Listing  (Continued) 


ill 

00  1^2  l-l.Nft 

1*^2  U?f  f  •  il*U?(  1 1  il<su«<4(  it 

PRlNT*,'*C0»'l>UTf  SU«  OF  ^Q\  OF  OlFfiftFNCFS  BtTwEEN  LOAOtO  ACTUAL  AN 
10  LOAOfO  PRfOICTCO  NOOAL  VfCtnRSs* 

c ••♦••##••#••«••••#••«• 

CA 

CA  READ  *^ASS-LCA0€0  NODAL  VECTOR^ 

CA 

CA#««*A 

00  ?S0  J-1*NC 
DO  ?50  I«1«NR 
UU<  f  f  JI’O.O 
READl^i?^! I  UUt 1 ,  i| 

261  F0RN4TI10».E16.^I 
260  PR INTA .UUI 1 • J) 

00  262  1-1«NR 
DO  262  J-lfNC 

262  UUt  I«  Jl  *001  I«  J}*uun  *  JI/IOO. 

DO  263  J«1*NC 

6UN1 1  JI-0.0 
DO  266  t-l«NR 

266  SUNK  J1*6UN1I  JKUUtI  f  JtAAe 
SUNl  (iKSORTISUNK  iM 
on  263  t-LiNR 

26  3  UUttf il«UU(r«i»/SU61(  il 

CA«A«A««A 

CA 

CA  CONFUTE  THE  SU"  OF  THE  SQUARES  OF  TmE  DIFFERENCES  BETWEEN  ACTUAL  AND  PRC- 
CA  OICTEO  LOADED  NODAL  VECTORS 

CA 

DO  263  J-l.NC 
OlFt JI-0«0 
DO  263  1-KNR 

263  OIF  I  JI-01F(  JKfUUt  I «  J 1  >U2f  I  •  i  I  I A  *2 
00  26Q  1*1, NC 

0IFKII«O«O 
260  OlFKlKOlFtll 

260  CONTINUE 

DO  261  l«l«NR 

261  OIFKII-SQRTtOIFlt  IM 

PRtMTA, -these  are  THE  VALUES  OF  THE  NODE  iHAPE  VECTORS:* 

00  266  1«KNR 

266  PRfNTA«(U2t (« it « i*L«NCt 
00  66  I«l«NR 
00  66  i«l«NO 
CC<1«  Jl«0«0 
0(7  66  R«1«N0 

66  CCt  !•  il*CC(  I  <  JMU3TP  1 1  *KlA|^lf  K«  Jl 
on  67  I«ltNR 

00  67  J-l «NC 
ANA! I, JI-O.O 
ANNfl, J|«o.O 
00  67  K«l«NO 

ANN|I  ,  Ji.amwu  t  i»*CCt  I  fKI  Au3Pf  K,  II 

67  ANAIT* Jl -ANNI I , it 
00  88  I«ltNR 

00  66  J-KNO 
00(I« JI-O.O 
on  88  K-ltNO 

86  00(1  •  ll•00t  I  «  il  «U3TPt  I  *1(1  AONECIIK  t If 

on  89  I«l,NR 
00  89  J-l,NC 
AAKtl, JI-0.0 
00  89  K-1,N0 

69  AAKt I • J I -A AK (1  , J  t  *rO ( I  .K |*U3P I Kf it 
on  72  r«|«NR 
on  72  l-KNC 

If(ANS<A-"<I,J»f,LC,lO.F-8» 

IF  1 A6S (A-A(|,J||.Lr.lO»E«Nf  AN»||,J|.0,0 
72  IF  (  A8S  t  AAPr  f  I  ,  j  M  .1  F  .  1  0«f '61  AAN  (  I  t  1 1  *0. 0 

CA 

CA  THIS  IS  THE  If AST-SOUARFS  SOLUTION  TO  THE  NASS  NATRIX  BEFORE  THE  NASS  HAS 
CA  been  AODEO  to  TMf  6,^  TF8«*. 

. . At 

Figure  C"l.  Progruni  ;iatrix  Listing  (Continued)* 


P«|NT*,“Tm|S  l*>  TH€  "*T8|X  knnt 
00  73  l-UNH 


73 

PR|NTP«**THIS  IS  TH(  MATRIX  AAK:** 

00  7A  l-l*hR 

7a  PR INTA, 1 AAK t I • JI • J-1 fNCI 

CAPPPPAP A* 

CP 

CP  SOLVE  THE  CCNFRALI7C0  FICENVALUf  PRORLfM  f OR  Tm€  “ASS-LOAOEO  NATURAL  FRE- 
CP  OuFnCICS  and  MOnC  shapes  using  the  rectangular  MOOAL  matrix  U3. 

CP 

CPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 

If  fKK.EO.ll  AMM|S,^l«AMAf^«MP0.lt4<i 
IFIKK.EQ.7I  AHM|S«SI«AmA|^«M»0«?A? 

IF  IKK.EQ.3I  AMMES«M«AMA(S.^|pO«P|S8 

17-9 

I  JQR-7 

lA-9 

ie-9 

N-9 

PRINTPt^THE  MOOIFIEO  MASS  MATRIX  AMR:* 

DO  67 

6  7  PRINTPt I AMM( I «  Jt • J-t  tNCt 

CALL  EIC7FfAAK«|  A,AMM,fA«N«|  J0P«AiEA7«Bf  TA7»/3,I7»t<Kp,ICR4,| 

CPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 

CP 

CP  THESE  ARE  THE  L E AS T-SflUARf S  SOLUTIONS  FOR  THE  NATURAL  FREQUENCIES  WITH  THE 
CP  The  AonEo  mass  in  rnt  s.s  positiqn  of  the  mass  matrix. 

CP 

CPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 

PR|NTP|*THf Sf  ARE  the  VALUES  Of  IHE  73  MATRIX:* 

no  7A6  l-UNR 

2kb  PR|NTP«(7  3<  I.JItJ-UNCl 

PRINTP. "THESE  ARE  THE  VALUES  OF  ALP:* 

00  IT  l-l,Nfl 

ALP( I l-CSORT I  ale A7f I  I /PETA7<  t  M 
HALPUI-CS0RT(ALPniPP7l 
31  PRINTP.ALPf n fWALPII  I 
CALL  VSRTAIWALP«LAI 

pRiNrp,*THfsf  ARE  The  values  Of  walp:* 

PR  |NTP,f  WALP<  M«  I-l  tNC  I 

on  76S  l-l^NC 

IFU.FQ.ll  AALP<9)-WALP<f  I 
tF<  I«CT« II  CO  TO  766 
CO  TO  765 

766  K-I-l 
AALPCKI-WALPIII 

765  PRlNTPfAlLPIKI 
00  19S  1-1  ,  R 
OQ  19R  J-1,NC 

IF f APS! AALPII |>RE AL <ALP( 111 t^LE .lO.E-ei  GO  TO  197 
IF  (AMS<  AALPI  I  I-REAL  T  ALPI  III  UCT*  lO,E-ai  CO  TO  190 

197  00  196  K-liNR 
UAfK,|  1-0.0 

196  UP|K.II«ll<>fK,II*KEALf73(K.flt 

198  PR  f NTP .UP  11  t  II 
on  iMo  i-i.NC 
SUMTI ll-C.O 

00  151  l-UNR 

151  SIJM3I  JI-SU-11  J»*UPU  .  IIPP? 

SUM3( JI-SORTI SUM3<  J|  I 
00  150  I-l.NR 

150  UPl I • Jl -UP (I . J I /SUMII ft 

PR  INTP  ,*CC*PUTf  THf  SUM  Qf  SOS  OF  niffCRINCES  HMVMN  10*011)  ACTUA 
IL  ANO  PRfOICTCO  LOAOEO  MODAL  VECTORS  USING  THE  GENERAll/tO  INVERSE 
7:* 

00  755  l-l.NC 
OIFFI JI-O.O 
00  755  T-I#NH 

755  OlFEIIt-PlFfl  JI«IUUlIt  It -UP  I  I »  II  |PP7 
DO  756  r-l .NC 

757  OtF7n  l-rif  f  I  I  I 

756  CONTINUE 
DO  767  l-l «NR 

767  01E7I f l-SCRTini E7C I  II 

.  Proqr.im  Matrix  Lij^tinq  ( ‘ '< -til  iriuod ) 


F’iqurc  C-1 


PB|NT*#"THf  €ICCNVKTOfiS 
on  29  1-ltNR 
PR|NT«,(U?( I • Jl« i-l iNCt 

PRINT***TH€  C£NEKAH/tO  INVERSE  EICENVCCIQRS  ARE:- 

00  79  I»liNR 

PRINTS* IUA( 1 • it • J«1 iNC t 

00  I-l iNC 

AALPI I t-*«lP( I  *  *100. 

UAMBOAI I |•MA^60A( 1 l«  L00» 
on  5?  L««ltNC 
REA0(9.S3}  N*t60AfLB> 

FORNAT(77X,Ee.3l 
pR|NTP*W**BOA|Lfll 
on  5^  1-1 .HC 

ERROR!  Il«(  (W*<eOA(  1  l-MA'^t^OAC  I  1 1  /W*<BOA|I  ||«100. 

00  ?3l  l-UNC 

ERROR  1  in*((H>1R0A<|t-AALPflt  l/N*tfiOAI  I  I  »•  100. 


C*  COMPUTE  the  nOT  PflOPuCT  OF  THE  ACTUAL  «ASS-LOAOfO  VERSUS  PREDICTED  MASS- 
CP  tOAOEO  modal  vectors  FOR  THE  SQUARE  MOOAt  MATRIX  IU7»  AnO  FQR  THE  RECTAN- 
CP  ULAR  MODAL  MATRIX  <U8I 

CP 


00  IPS  I*ltNR 
00  IPS  J-l,NC 
USd,  JI-0.0 
U6f I«  ]l«0.0 

US!  t • JI«US(  !• JI*U2<  J« I > 

IPS  U6f  t  t  JI-UM  t  •  Jl  «UU!  J«  r  I 
QQ  IP6  t«l«NR 
00  IP6  J-1*NC 
U7! f , JI-0.0 
UM  !  I • Jl •O.O 
00  R«l*NC 

U7! I  •  ll«U7t ( I J1 «UU! I »RI PUS!K«  H 
IP6  U8!I«JI-ue!T«JdUU!r  «t<|PU6IK,  if 
on  IP7  I«1,NR 

IP?  PPtNT!f*.lP«P»  !U7((,f»,J«l,NCI 
UU  IPP'*  1*1, NR 
IPP9  PPfNTCSi  tPPMI  <UM!I, 

IPP8  F0PMATt7X,9!2X,F 7, 1 1  I 

PRTNTP^-'THE  ElGENVALUCS  ARE:** 

PRiNTPfMTHF  COMPARISON  OF  E  tCENVALUES  :  ** 

PR|NTP»**  UNIOADEO  PREOtCTCC  ACTUAL 

IRE  PPEOICTEO  percent  SQUARE** 

pRiNTPf**  experimental  -ASS-LOAOFO  MaSS-LOAOED 

it  mass-loaoeo  error  root* 

PRINTP*"  DATA  USING  EXPERfMENTAl 

I  USING  PREOICTFO  of** 

PRINTP,**  SQUARE  MODAL  DATA 

1*5  PSEUOOINVERSE  VERSUS  OIFF*S** 

PRINT*,**  matrix 

10  !9RnwS,8C0LSI  ACTflAt  SO*D- 


PERCENT 


EXPERIMENTAL  PREDICTED 


PRINTP,**  SQUARE  MODAL  DATA  VERSUS  OIFF 

1*5  PSEUOOINVERSE  VERSUS  OIFF*S** 

PRINT*, •  matrix  actual  SO* 

10  !9RnwS,8C0LSI  ACTflAL  SO*D- 

00  29  I-l,NC 

OMEGAS! If •OMEGAS! I  1*100. 

PRINT!  ft, SH  0“EGAS!I  l,N*“POA|  1 1  •WMBOAf  1)  ,  E  RPQR  I  I  I ,  D I  F  U  I  I,  AALPf  I  I  • 
lERRORl ! 1  I ,0IF2 ! 1  ! 

FORMAT!PX,F7.3,7K,F7,T,7X,F7.3,ftX,F6,e,PX,FS.7,SX,F7, J,7X,Ffe.^,PX, 

IF  S  #  ?  1 

TFILLL.EQ.il  CO  TO  2b9 
IFfLLL.EO.il  GO  TO  ?ft9 


C*  IF  LLL  EQUALS  ONE,  THIS  5TFP  WILL  PC  IGNORED  AND  Thf  PPOGRAM  MILL  CO  TO 
CR  STATEMENT  ?ft9.  HOMEVER,  I*  LLL  FQUALS  7,  The  MASS-LnAOED  "DOAL  VECTORS 
C*  lACTUALI  MILL  BE  READ  IN  ANO  TME|R  mqDE  SHAPES  PLOTTED, 

CP 


2S8  00  768  J«t,NC 

on  ?#>8  I«1,NR 
REAniS,7ft7l  U7fl,il 

767  F0RMAT!nx,F:0.7l 

768  U7II,  II -t^?!  I,  if  7  100. 


Figure  C-1.  Program  Matrix  Listing  (Continued)  . 
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00  1^0  J-l«NC 
J»-0,0 

00  141  1«l*NR 

141  SU4?|  JU5U**?1  !»•♦? 

SU42I J|-SORTf$UH2l J) I 
00  140  t«l.NR 

140  U2|I«n-U2II,JI/^UP<2(  Jl 

269  CONTINUE 

on  211  LL-l»NR 
CALL  RCNPLIll 

CALL  TITL3ni««00€  tCONP  f  CUft4  T  t  ON  i*» ,  100  *  4 . 0 . 4 . 0 1 

AXC^^nt«TV4£TA%«tlOO*«ri'**lOO***«AOl4L4**lOO*l6.  tlS.  «  16«  t 
CALL  VUA4SI-20.*70,.S0.I 

call  CRAE30t-20, *l0.»^n.*-l«,10**60«*-0.6*l0. *^0. \ 

LK-LL 

call  NFMPfNtll 
00  21?  KL-lt2 
IFlKL.PO.il  CO  TO  ?04 
IFIKL.CT.II  CO  TO  ?0l 

204  00  206  1-1.3 

?06  CALL  CU0V3P1K |K| ,YtKl ,70IKI •3.1 1 
00  ?07  J-1,3 
K- J 

VV1ll-K|Kt 

?2I11-?01KI 

XY|2I-X(K«3I 

YY12I-Y(K.3I 

??f?l-?0(K*3l 

XX|3I-Y|K*6I 

YYI3I«Y|K«61 

2?|3l-?ntK*6l 

PR|NT*.YY|ll.YY11|,/;ill.KX12l,YY1?).Z?|2l.)(Xt3I.YV13l.?2l3l 
207  CALL  CUPViniYK.YY, 21,3.11 
IFIKL.PO.ll  GO  TO  203 
201  on  203  L-l.NR 
YYYILI-XfLl 
YYYai-YILI 

IFILLL.PO.n  CO  TO  199 
IFILLL.CQ.2I  CO  TO  199 
IF1LLL.CQ.3I  2/21LI-^0<LUU41L,LLI43. 

CO  TO  203 

199  /?21L>*/Oa»^U2U.LLIA3. 

203  PRlNTA.XIXlLlfYYYlLl f Z221LT 

205  00  209  1-1,3 
K-3AtI-ll*l 

209  CALL  CURV301XXX1KI,YYY|Kl,r221KI,3,H 
00  210  J-1,3 

K- J 

XAIll-XYVfKf 

Yifll-YYYIKI 

ZA11l-T771Kt 

X4I2I-XXX1K43I 

YA12I-YVY1K-3I 

2AI7I-Z721K43I 

Xi13l-YXX(K«6I 

YAI  3l-YYY|K4f,| 

ZAm-7rriK»6i 

PRINT«.YAUI,YACl|,2AU»*XA12l,Y«f2l,2AI21,XA13).YA|3t,2Ani 

210  CALL  CURV30IXA,VA,7*,3,l| 

IFIKL.CT.II  CO  TO  216 

20«  CALL  NnMIOF 

CALL  NFMP£N|1» 

IF  1KL.FQ.ll  CG  TO  217 

216  CALL  NFWPFNUI 

217  CONTINUF 
CALL  0A3H 

212  CONTINUF 

211  CALL  fNOPLIU 
6  CONTINUE 

CALL  OONFPL 
215  STOP 
CNO 

'  FNO  OF  trcn®n 


Figure  C-1,  Procji  cim  >'atrix  Listing  (Continued), 
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-19.  l?.^ 

17. 

19.  ?4,37S  10.^ 

-19.  6. 

-9.9  6.  l?.9 

19.  6.  10. 

-19.  0.  1?.9 
-9.9  0.  1?.9 

• 

1  19.  0.  10.9  1 

1 

1 

i<,<i.44r  7 

.1969699 

9 

1 

Z  90 

2 

141.781  7 

•2607813 

9 

7 

1  103 

0 

706.  I*** 

,7970299 

9 

1 

9  103 

0 

706.641 

.7313899 

9 

9 

2  119 

0 

714. 010 

.6190009 

9 

2 

9  tl9 

0 

714.717 

. 7679099 

9 

9 

9  139 

1 

279.993 

,7939196 

9 

9 

2  229 

1 

991.669 

.7763902 

9 

2 

3  29? 

1 

484.447 

.7067292 

9 

3 

1  1 

160. 

.OOOf *00 

.OOOf *00 

30.9 

.0 

.0 

1  2 

7.69 

.OOOf *00 

.OOOf *00 

72.6 

•  0 

.0 

1  3 

-29.6 

.ooot»oo 

.OOOf  *00 

99.9 

•  0 

•  0 

1  9 

93.3 

.OOOf *00 

.OOOf  *00 

63.  1 

•  0 

.0 

1  9 

-99.2 

.OOOf *00 

.OOOf  *00 

69.9 

•  0 

.0 

1  6 

-31.2 

.OOOf *00 

•OOOf^OO 

22.  7 

•  0 

•  0 

1  7 

10.9 

, OOOf *00 

•  OOOf  *00 

106.9 

.0 

.0 

1  0 

9.10 

.OOOf *00 

.OOOf  *00 

99.1 

.0 

.0 

1  9 

-11.3 

.ooot^oo 

.OOOf  *00 

90.  7 

.0 

.0 

2  1 

132. 

.OOOf *00 

.000E*00 

121.9 

.0 

•  0 

2  2 

2.91 

.OOOf  *00 

•  OOOF  *00 

3 

•  0 

•  0 

2  3 

-27.6 

.OOOf  *00 

.OOOF  *00 

169.7 

.0 

•  0 

2  9 

127. 

,ooor*oo 

•OOOf *00 

109.  1 

.0 

•  0 

2  9 

-99.9 

.OOOf *00 

•  OOOf  *00 

122.9 

•  0 

•  0 

2  6 

-23.9 

.OOOf *00 

.OOOf  *00 

28.  1 

•  0 

*0 

2  7 

9.91 

.OOOf *00 

.000E*00 

172.1 

.0 

•  0 

2  0 

0.39 

.ooor^no 

.0006  *00 

9.  3 

.0 

•  0 

2  9 

-19.9 

.OOOf *00 

.0006*00 

139.2 

.0 

•  0 

1  1 

00. e 

.OOOf *00 

.OOOf  *00 

177.1 

•  0 

•  0 

1  2 

-19,6 

.OOOf *00 

.OOOE*CO 

122.1 

.0 

.0 

1  3 

193. 

.OOOf *00 

.0006  *00 

171.8 

•  0 

.0 

1  9 

790. 

.OOOf *00 

.OOOf *00 

191.6 

•  0 

*0 

1  9 

202. 

.OOOf *00 

.OOOF  *00 

190.9 

*0 

•  0 

t  6 

-780. 

.OOOf *00 

.OOOf *00 

139.8 

*0 

*0 

1  7 

-192. 

.OOOf *00 

.OOOf  *00 

117.1 

*0 

•  0 

1  9 

-29.1 

.OOOf* 00 

.0006*00 

170.1 

.0 

•  0 

1  9 

193. 

.OOOf *00 

♦  OOOf  *00 

120.6 

.0 

•  0 

9  1 

90.9 

.OOOf *00 

.0006*00 

99.2 

•  0 

*0 

9  2 

-11.0 

.OOOf *00 

•  OOOf  *00 

111.8 

*0 

.0 

4  1 

148. 

.OOOf *00 

.OOOf *00 

106.9 

•  0 

•  0 

9  9 

737. 

.OOOf *00 

♦0O0F*OO 

139.'. 

.0 

*0 

9  9 

226. 

.OOOf *00 

.oobf *00 

131*9 

•  0 

•  0 

9  6 

-279. 

.OOOf *00 

.0006*00 

126.  7 

•  0 

.0 

9  7 

-121. 

.OOOf *00 

.0006*00 

122.6 

•  0 

•  0 

9  9 

-29.9 

.OOOf *00 

•OOOf *00 

119.2 

*0 

.0 

9  9 

190. 

.OOOf *00 

.0O0E*00 

113.9 

.0 

.0 

2  1 

269. 

.OOOf *00 

.OOOf *00 

69.  9 

*0 

.0 

2  2 

-19.0 

.OOOf *00 

.000E*00 

190.2 

.0 

•  0 

2  3 

192. 

.OOOf *00 

.OOOF*00 

191.9 

.0 

•  0 

2  9 

77,2 

.oooe*oo 

.0OOE*OO 

96.9 

.0 

•  0 

2  9 

-131. 

.OOOE*00 

•OOOf *00 

199.9 

.0 

.0 

2  6 

-902. 

.00OE*00 

.0006*00 

190.0 

.0 

•  0 

2  7 

-96.9 

.OOOf *00 

.0006  *00 

199.  7 

.0 

.0 

2  6 

-I'^.O 

.000f*00 

.OOOF  *00 

199.  7 

.0 

.0 

2  9 

1 12. 

.ooor*oo 

.000£*00 

192.6 

.0 

.0 

9  1 

396. 

.OOOF *00 

.000F*00 

96.9 

.0 

•  0 

9  2 

0.91 

.0007*00 

•  OOOF  *00 

29.0 

.0 

•  0 

9  3 

07.7 

.OOOf *00 

.OOOf  *00 

161.9 

.0 

•  0 

9  9 

33.6 

.OOOf  #00 

•OOC! *00 

96.  7 

.0 

.0 

9  9 

-ICl  . 

.OOOf .00 

,000f  *00 

190.6 

•  0 

.0 

9  6 

-917. 

..1001  ♦<'0 

•  OOQF 

•  17.9 

•  0 

.0 

9  7 

61.6 

,ooo*^*oo 

.OOOF  *00 

1  3.9 

•  0 

•  0 

9  9 

-2*’. 6 

.loor  *00 

•  OOOf  *1)0 

191. n 

•  0 

•  0 

9  9 

03.  1 

.000F*01 

•  OOOf  *00 

160.  3 

•  0 

.0 

9  1 

923. 

.OOOf *10 

•  OOOf  *00 

22.  3 

•  0 

.0 

9  2 

-11.9 

.OOOF  *00 

.OOOf  *00 

121.9 

.0 

.0 

9  3 

129. 

.OOOf  *00 

•  OOOf  *00 

1  79.9 

•  0 

.0 

9  9 

00.0 

.006F.OO 

.OOOf  *00 

21.1 

•  0 

.0 

Figure  C-1. 
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4 

3 

-4C.8 

.OOOr  *00 

.ooot  ' 

64,9 

•  0 

.0 

4 

6 

737. 

*00 

.OOOf  ♦«  .' 

i  89,8 

.0 

.0 

4 

7 

87*8 

.OOOf *00 

•  OOOf  *0f^ 

M.4 

.0 

•  0 

4 

e 

•^4«0 

.000€400 

.OOOf  *00 

173.8 

.0 

•  0 

4 

9 

107. 

.0007*00 

.0006*00 

174.  7 

.0 

.0 

2 

1 

-774. 

.OOOf  *00 

,0006*00 

1  77.  7 

.0 

•  0 

2 

2 

-77.7 

♦  OOOf  *00 

•  OOOF  *00 

88.0 

.0 

•  0 

2 

3 

701. 

.OOOF  *00 

.0006*00 

36.3 

.0 

.0 

2 

4 

33.7 

.000€*00 

•0006*00 

116.8 

•  0 

•  0 

2 

3 

88.0 

•  0006  *00 

.0006  *00 

133.6 

.0 

•  0 

2 

6 

-30.8 

•0006*00 

.0006*00 

7.9 

.0 

.0 

2 

7 

-89.3 

•0006*00 

.0006*00 

94,7 

•  0 

.0 

2 

6 

-18.0 

•0006*00 

.0006*00 

6.3 

.0 

•  0 

2 

9 

16.0 

•0006*00 

.0006*00 

14.6 

•  0 

.0 

3 

1 

979. 

•  0006  *00 

.0006*00 

L  1.  1 

.0 

•  0 

3 
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Figure  C-1.  Program  Matrix  Listing  (Continued). 
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figure  C-1.  Program  Matrix  I.isting  (Continued). 
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Figure  C-1.  Program  Matrix  Listing  (Continued). 
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9 

7,94 

.0007*00 

•  OOOF 

*00 

121.4 

•  0 

.0 

7 

6 

-133. 

.0007*00 

•  OOOE 

•  00 

144.4 

.0 

•  0 

7 

7 

-21.6 

.0007*00 

•  OOOf 

♦00 

198.4 

.0 

•  0 

7 

a 

41,9 

•0007*00 

•  OOOf 

•  00 

191.0 

•  0 

.0 

7 

9 

99.4 

.0007*00 

•  OOOE 

♦  00 

136.2 

.0 

•  0 

9 

1 

244. 

.0007*00 

•  OOOE 

•  00 

11.8 

•  0 

.0 

9 

2 

-44,0 

. 0007*00 

♦  OOUf 

♦  no 

7.4 

•  0 

.0 

9 

3 

-94.8 

•  OOOf  *00 

•  OOOF 

•  CO 

4.  7 

•  0 

•  0 

9 

4 

-103. 

•  OOOf  *00 

•  OOOf 

♦  00 

144.8 

.0 

.0 

9 

9 

-19,4 

.0007*00 

•  OOOE 

♦  00 

147,9 

.0 

•  0 

9 

6 

144. 

•  OOOF  *00 

•  UOOf 

♦  00 

173,4 

.0 

•  0 

9 

7 

-84.9 

.0007  *00 

•  OOOF 

♦  00 

1  74.2 

.0 

.0 

9 

a 

-72.3 

.0007*00 

•  OOOF 

•  00 

.8 

.0 

.0 

9 

9 

-94.2 

.0007  *00 

•  OUOf 

♦  00 

•  1 

•  0 

•  0 

6 

i 

977. 

•  OOOF  *00 

.OOOF 

♦  00 

91.7 

.0 

.0 

Figure  C-1.  F^ro  rrani  .'iatrix  Listing  (Continued)  . 
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6 

2 

28. 

•  OCOf.  *00 

.coot •  0 

176.9 

.0 

,0 

6 

3 

-33.8 

.OOCE  *00 

•  OOff  *‘-  i 

106,  3 

•  0 

.0 

6 

4 

-42.2 

•  OOOf  *00 

.OOOf  *00 

109,  1 

•  0 

.0 

6 

5 

4.72 

.OOOf *00 

•  OOOf  *00 

29.4 

.0 

.0 

6 

6 

-41.1 

.OOOf  *00 

.OOOf  *00 

29.  1 

.0 

.0 

6 

7 

-188. 

,000€*00 

•  OOOf  *00 

178.1 

•  0 

.0 

6 

8 

-24.2 

,ooor  *00 

•  OOOf  *00 

42.6 

•  0 

.0 

6 

9 

-74,0 

.OOOf  *00 

•OOOf *00 

72.2 

•  0 

•  0 

7 

1 

432. 

.OOOf  *00 

.OOOf  *00 

6.6 

•  0 

.0 

7 

2 

-274, 

•  OOOf  *00 

.OOOf  *00 

28,1 

.0 

.0 

7 

3 

-464, 

.OOOf  *00 

.OOOF  *00 

,  7 

.0 

7 

4 

448. 

.OOOf *00 

.OOOf *00 

14.4 

.0 

.0 

7 

8 

-14,4 

.00  E  *00 

.OOOf  *00 

14,9 

.0 

.0 

7 

6 

-no. 

.OOOf *00 

•  OOOF  *00 

14.3 

•  0 

•  0 

7 

7 

442. 

.OOOf  *0 

.OOOf  *00 

27.3 

.0 

.0 

7 

8 

-223. 

.000*^  *00 

.ooor  *00 

l.l 

.0 

.0 

7 

9 

420. 

•  OOOf  *00 

•  OOOF  *00 

167. 3 

•  0 

•  0 

3 

1 

344, 

.OOOf  *00 

.OOOf *00 

1  12,0 

•  0 

•  0 

3 

2 

7.72 

.OOOf  *00 

.000E*00 

18.3 

•  0 

•  0 

3 

3 

43.8 

.OOOF  *00 

.OOOf  *00 

99.3 

.0 

•  0 

3 

4 

-88,3 

.OOOF  *00 

.OOOf *00 

162.1 

.0 

•  0 

3 

S 

-3.97 

•  OOOf  *00 

.OOOF  *00 

34,4 

.0 

.0 

3 

6 

14.4 

.OOOf  *00 

.OOOf  *00 

34.9 

.0 

•  0 

3 

7 

-44.0 

•  OOOf  *00 

.OOOf  *00 

21.4 

•  0 

.0 

3 

8 

20.6 

.OOOf  *00 

.OCOf  *00 

40.0 

•  0 

.0 

3 

<9 

24.1 

.OOOf  *00 

•OOOf *00 

73.2 

•  0 

•  0 

3 

84 

2 

168.284 

2.8288817 

4 

3 

90 

0 

180.711 

.6866900 

4 

4 

108 

0 

204.393 

.7893097 

4 

4 

6 

110 

4 

220.972 

4.8140106 

4 

6 

7 

119 

1 

238,822 

1.0497084 

4 

7 

5 

139 

0 

279.744 

.4410013 

4 

4 

6 

144 

1 

1.0753403 

4 

6 

7 

150 

0 

300.294 

♦4929404 

4 

7 

3 

142 

4 

2.8444467 

4 

3 

m 

ewo  Of 

«fC0«0 

l-fNOt 

Figure  C-1. 


Program  Matrix  Listing  (Concluded), 


101 


^66.C76 


-2.M9 

-  .  2  C  2 

1 ,  f  ^ 

1.4)7 

3‘6 

-  .  1  2 

.  1 

.  .'8  * 

-4, 310 

-  .  L'  6  7 

,  6  "  4 

-.12? 

2.862 

-•Ul4 

- .  -  n 

2.818 

793 

“  .  7  2  8 

.  244 

.282 

.407 

IHl 

•  S04 

-.028 

732 

-,C  39 

.43? 

—  •  i>  6  6 

.  162 

.00  8 

-.099 

UU  ? 

-.861 

.0  26 

.  120 

-.086 

52,106 

2.66  8 

22  .903 

-4,431 

2.666 

.  173 

1,438 

-.3  00 

22.903 

1.408 

12.884 

^2.98  ? 

-4.431 

-.300 

-2.96? 

4 . 2h2 

.909 

.090 

1 .108 

.  19Q 

-6.920 

-.263 

-3. 789 

-1.661 

10. 783 

.580 

4.842 

-1.821 

-2.448 

-.117 

-1 .  IC  ? 

.  173 

4.353 

.265 

2.625 

-.20  7 

The  EICENVALUPS  ACE: 

The  C0-PA9ISQN  OF  ClGtNVALUES: 

UNLOADFO  PCfOlCTEn  ACTUAL 

EXPERI1ENT4L  HASS-LUAOED  '^ASS-L0A0£D 

0*T*  USING  firPf  e  1  Hf  NT  At 

-  -  -  SOUAke  SCIUL-.  -  OATA 

•MATRIX 


lfc9.A37 

206.396 
206. 6A3 
—  .Z3a.039 
238.212 
279.993 
A51.669 
989.A87 
I. 


l  ?6  .  'Jti  A 
197.687 
206. a99 

.233.101 _ 

238. lAA 
278.672 
A^.0. 266 
^83. 779 


1 39.97S 

1  7Q, fl32 
20A.2^^ 
20A, ftSS 
238.630 
238,660 

2  60 . 6<»6 
AS3.66a 
466.076 


*Q393526  32  761Q<, 
-.178363006826 
•  2073iq?0'J8766 
-.  37639248  7;q^ { 


-.  16701314«6U6 
♦09739424924918 
.03701173771741 
-.1063298400984 

.624l90O77S4.fl6 
-.048704^6  746  189 


-.3588833139861 

1. 


r. 8188871713678 
.0246809«0«6104 

.I4149l2il»^3^  , 

•  1  12  1  340266^6? 


-.37646369 

.2970937382'^‘'7 

-.0Q24C966249?2t6 

.065943013791/7 

1. 

•  1533887643867 


I-.244C58  34  78  7/2 

j-. 1 18260 1321093 
[-.000M422  366S62R7 


*■  i  i 

.  l4t. 

o  6  9 

J6/ 

f  u 

-.141 

/82 

.231 

462 

.003 

2  1  i 

-1.097 

2  1  f 

361 

u  6  0 

-.1  16 

4  v9 

-.001 

•iKji 

-9.920 

0  90 

-.261 

10a 

- i. 789 

190 

-  1 . 6M  1 

dr.  3 

-.d68 

868 

7.531 

Olo 

-.173 

184 

•  868 

al2 

-1.941 

.f  -i  i 

-.281 

-  .  -  1 

-.022 

-  .  ^  j 

-  .  04u 

.12? 

•  L06 

-.116 

-.038 

-1.12«, 

ia8 

.Id? 

-.081 

-.dVi 

.026 

.049 

.001 

lo. ?a J 

-2 . 44d 

.  8  8u 

-.117 

4  .  d  42 

-1.107 

-1.821 

.173 

J  iO 

-.  Id4 

1  73 

.865 

2.697 

-.486 

-.4  66 

188 

.6l4 

-.270 

Il4 

•  c  *.?> 

-.u2l 
-  •  2  Oo 
.182 
.66d 
.03  2 
-.018 
.  008 
4.383 
.268 
2.628 
-.207 
.812 
-1.94j 
.614 
-.270 

•  904 


percint 

SOUARt 

PREDICTED 

PFRCENT 

SQUARE 

tKRUR 

ROOT 

‘lASS-LDAOeO 

Error 

KOJ  T 

PRfOKTFO 

OF 

U$  I  NO 

PREOICTEO 

OF 

versus 

Olf f  ’S 

PSfcUDO INVERSE 

VERSUS 

OIFF  •$ 

actual 

so»o 

1 9R0wS . 8C0L  S  > 

actual 

80‘D 

-13.92 

.39 

169.437 

-21.08 

2.59 

1.82 

.88 

ldl.261 

-.79 

2.65 

3.21 

.39 

206. 396 

-1.08 

1.44 

-.90 

♦  42 

206.643 

-.97 

1.49 

-  2.32 

3.fi7 

23B.039 

.25 

3.00 

.22 

1.31 

238.212 

.19 

3.27 

•  63 

.84 

279.993 

.23 

1.26 

•  76 

0,66 

481 . 669 

.  44 

8.71 

-3.80 

.72 

0.000 

lOO.CO 

1.  79 

Figure  C-2 


Edited  Output  From  Program  Matrix 


S.'.niplo  output 


ifir.' 


'roin  ^pvtpiiii.-^. 


hne'<  an  <-r  I  S'  j 


V  in  a  If  I  S  *  I  •  v*i- 
y  in  A  »  I .  i  .  >:t 

!N  A-s  .  3-^  ‘iS  I  rs 


VI FMPO  !NT 


xvu>-2.coce*Li 
YVU-  7.000?  *0  1 
ZVU*  5.000( «01 

IN  AfiS.  3-0  UNITS 


C8APM  Sf T-UP  (  CfiAf  30  I 


ilRlCIH  _  . 


X30nBlClN--2.000£*0l 
Yinno iciN«-i .oocc  *00 
73Dl)R!CIN— S.OOOE-Ol 

STEP  S176  .  _ 


X10STP*  l,0uyE*Cl 
Y3nSTP*  I, 0006*01 
Z3DSTP-  l.Cu06«01 

*^1X1  •u- 


vin-Av  0006*01 
Yin'*4x-  s. 0006*01 
730’*AJ(-  S.00Qt*0l 


LOCATION  OF  CU8>J6NT  PHYSICAL  OfilClN 

X-  .50  y-  1.12  iNChes 

ppfi’*  tnwffi  LfFT  cnPNfp  nr  pace 


Figure  C-2.  EcJitod  Output  From  Program  Matrix 
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401.47a 
434.504 
456. 5?4 

^.40^ 
.137 
1  •40Q 
1 . 0  a 
-.005 
1.713 
-.7fl7 
.  1C<J 

17.  3<»3 

-•64  3 
-1.041 
.45  3 

-.450 

.  - 

7.7  16 
-.010 
1.779 


-.141 

.  14  7 

1...U 

—  .67*, 

1.474 

.  7  n  .1 

.  4  79 

-.618 

,  0  Ih 

-,u7l 

-  •  0  3  6 

-•VJ64 

7  7h 

.074 

-.o33 

.1..75 

•  05ti 

-.447 

.  567 

-.  08  7 

.08  3 

•  52  4 

•  ul2 

-.163 

-  .  (M  ‘‘ 

-.  J4.) 

1.605 

.  i  5tj 

.  147 

.494 

.  049 

-.160 

-.041 

.055 

-.1  16 

•  7  82 

.  c  8  7 

.  j  39 

-.037 

•  00  8 

-.178 

-.141 

-7.930 

.861 

-7.  J19 

.4  71 

-.717 

-.037 

.044 

-.157 

»OHI 

-.374 

-.1  39 

.  Ou  7 

.  L  18 

-.2  34 

.074 

-.<'54 

.079 

-.115 

-.  164 

.005 

.004 

.001 

*•001 

“  .  1 6  9 

.657 

-.047 

.8  79 

.454 

•  118 

370 

-.f  4  1 

-1.041 

.44  1 

-.450 

-6.777 

7.736 

-.010 

1.  779 

.137 

.7  72 

-.013 

-.0  75 

.  1  75 

-.C4l 

•  065 

-.098 

.77  7 

1  .611 

1.6  19 

-•067 

358 

-.4  73 

.073 

.967 

-.01  3 

1.639 

7.867 

•  0  7a 

-1.879 

-.349 

-.040 

1.594 

-.075 

-.067 

.076 

.794 

.943 

4O0 

-.137 

.077 

•  1  75 

- -.35A 

*1.879 

.943 

9.367 

-'1.835 

-.307 

-1.747 

-.041 

-.473 

-•  149 

-.400 

-1.835 

1.757 

.  159 

-.047 

.065 

•  073 

-.040 

-.132 

-.307 

.  159 

.077 

-.048 

-.098 

.967 

1 . 594 

.072 

-1.747 

-.047 

-.048 

1.214 

5  AflF  : 

—  UliiLQAQtO  __Pft£CICTtQ 
f)fPpffr"fNTAL  ••A5S-UnAOfO 
data  U5tNC 

50U4«E  4Q0AL 
TtATtflX 

137.775  137,085 


137.775 
—  207. oa& . 
739.187 
760.553 
794.845 
377.881 
394,774 

- 414..1A4- 

453.688 


ACTUAL 
••A55-LOAOfO 
eXP£PI IfcNTAL 
OATA 

149.877 


--Ul.iOO- - 707.598 

777,943  738.415 


760, 788 
786,7  78 
377.150 
393,434 

- 413.4^4 

453.491 


776.147 
301.519 
373.745 
401 ,5  78 

- 434.504 

456.574 


.008457103553768 
-.1817037949791 
-.I973675p^7q'»3 
-.04608525837178 
.22681980  344,47 
.1959742317667 
.01445647107975 
-.1195120697767 
.03897287641419 
-.0513577429707 
.2057643600871 
1  . 

.5553ll94lOf«81 
-.7756695747378 
-• 17914404  7641 1 
-.07121 9  166/4-148 

.230377*717447 
-.324163584 7 5»1 
.09765044P  7  '9'’ 
-.70C52173675!! 
.7396044  777/47 
-.44  1  37  60  8‘<^4  7? 


P£fiC£NT 

SQUAR£ 

P&ED1CT£0 

PEftCf NT 

SOUARi 

EfiPOP 

P90T 

•^ASS-LOAOED 

EPROB 

ROD  r 

»R£D1CT?0 

OF 

USING 

PPEOICTEO 

OF 

VERSUS 

OIFF*S 

PSEUDO  INVF.PSF 

VF  RSUS 

OIFF  »’ 

actual 

SQ*0 

19RQ4S.aC0LS I 

actual 

SQ*0 

8.53 

.45 

137.775 

8.41 

.43 

-  10.56 

.  .44 

-  207.085  ..  . 

-  -7.21 

2.57 

4.39 

3.96 

239.  187 

*.  32 

4.28 

5.56 

1.85 

280. 553 

-1.60 

3.  04 

5.05 

1.67 

794.848 

7.21 

2.53 

•  1.71 

♦  77 

377.881 

-1.43 

•  89 

7.03 

.17 

394.724 

1.83 

•  35 

-  4.74 

—  *  16  - 

-  414.164 

4.68 

.26 

.66 

4,08 

.001 

100.00 

7.99 

Figure  C-2.  Edited  Output  From  Program  Matrix  (Continued) 
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?7Q,  7b*> 
1  7S 


-.628 

2  .  "  6  ^ 

.  ? )  7 

1  .  6  6 

•  2  4C 

sno 

\  •^*'2 

•  •• 

“  .  20  0 

.041 

-.99#. 

-.749 

1.814 

.  7^> 

-.117 

-  .  2  7 

-.^os 

-.131 

,  6  nf» 

.  7  ■)  » 

1.0)2 

.  1  /s 

-.2  7« 

.0*3 

-.140 

-.09) 

.0*4 

-.107 

1.167 

.02  1 

-1  .  F4N  \ 

-  ,  f'  ,4 

-.119 

1  •  Ol;M 

.333 

1  .  9  4 

-  .  L-  1  O 

.  )  7  0 

.  L  80 

-.824 

-.  397 

720 

.  'SI 

1  76 

..  H2 

-  4  #.  4') 

-.611 

.192 

1.91J 

1.142 

.2  16 

-1.220 

7. *03 

1  ■»o 

-1  .  '  M 

-1  .OSC' 

.08  1 

-  1 .241 

-.lie 

.702 

1.072 

-.629 

-.Oil 

-.887 

-1.371 

1 .072 

2.1)7 

-.439 

-.046 

-  l.OQ) 

-1.990 

-.620 

-.  4  )0 

2.089 

-.09b 

.387 

.081 

-.01 1 

-,U46 

-.rtS6 

.0  36 

•  004 

-1.241 

-.487 

-l.C'O) 

.  187 

•  004 

2.96) 

.877 

-.292 

-1  .  160 

.086 

.i»  1  7 

-.9  »2 

-.064 

•  645 

.970 

-.997 

•  009 

-.921 

.♦'9  .1*?  .0?.’ 

EICENVJLUf^  ARf: 

CP"PJ»!SnM  nf  f  1  rF  MVAL  uF-i : 

t  .CIO 

-.069 

-1.9)9 

,  I'J-J 

. 

.7^1 
,0U 
.  /  J7 
.  714, 
»2'>b 
.  4.^2 
.87? 
.2^72 
.  IS*^ 
•  U  9b 
.017 

.90) 

.184 


.<7b" 

.  iJ  0 1 

-.  J21 
.  J2H 

-.0^4, 
-.  7b4, 

.179 

.00'? 

.747 
-.064 
.649 
.  >?  74? 

-.997 
.00*i 
-.  ?21 
-.  184 
.621 
.241 


1.04) 

•  14»6 
.261 

?‘? 

-1.6  ?0 
.2  >9 
.191 

1. 944 
.499 

•  142 

•  072 
l.OlO 
-.06? 

-1.939 

.‘?94 

.241 

1.84*? 


UNLOAOEn 
EVPEff IMfNTAL 
.  QkTA 

167.266 
180.424 
207.298 
228.945 
238.930 
27<?,694 
2^<?.232 
301.621 
328, 9?3 

•07796039991111 

•  96<*17f  oco*,;4? 
-•4l64ll86fl92l‘- 
• 324C401 102942 

1. 

324<»4l96940g4 

•  199193<?140799 
•.1024289017242 
.03649737474271 
.3779160106226 
1. 

-.9996494967106 

•  464  7784  1  704  3 

.0l894646Bni^« 

3441949906034 
.  197241  2U2'’4  1 
14784  76192  76 
•.2729614399000 
•  96  71664094602 
-.8309980742906 
1. 

•  787999900419 
11474)0736796 
.1441790022 ?7’ 


POf niCTED 
•499-L^AOf 0 
USIVC 

90U4PE  ■00*1. 
■ATQ  I  V 
109.091 
179,69b 
199.339 
227,745 
237.442 
262.269 
290,98) 
101.309 
129, 30-? 


*f  TUAL 

PEOCE  M 

SOUAPf 

PPEOICTEO 

PERCENT 

SOUAftE 

A99-LOAnf 0 

1  ftp  09 

ppn  f 

■A99-LI)AOFn 

E  PQOft 

ROOT 

XPf  «?  I  “ENTAL 

PPEOICFEO 

OE 

USING 

PPfOICTf  0 

CE 

OAf  A 

VEft9U9 

OlfE  •$ 

PSFUDOiswteSf 

Vf  BSUS 

01  EE  *9 

ACTUAL 

$0»D 

l9ft0wS»8C0LS» 

ACTUAL 

so*o 

1  64.. >89 

37.99 

J  .87 

167.266 

,61 

1,93 

180.711 

.9H 

2.14 

180.424 

•  l6 

2,21 

2'J8,  365 

4.34 

1.26 

207,298 

.91 

2.  16 

220. 07^ 

-3.06 

1,91 

228,943 

-3.43 

1,60 

2 )8.822 

.98 

1  .64 

238,930 
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Figure  C-2.  Edited  Output  From  Program  Matrix  (Concluded) 
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Figure  C-2.  Sample  NIASTR.^”^  Dock  (Continued)  . 
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Figure  C-2.  Sample  NASTRAN  Deck  (Concluded). 
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Fiqure  C-3.  Data  Acquisition  Figure  C-3.  (Concluded). 

Program  for  HP5451B 
Fourier  .\nalyser 


TEST  I[€HTIF1>T10N  ISBQF  tAVE 

Ci'*l1A0HEnT  X,V,Z  tJAICIN  ORIEHTATION  CODE 

f’f?  1  OHO*  .e»»»  18  3* 


Sample'  Y-9  Modal  Analyf^is  Set  Up  with  A.)  General  Information 
Identification,  and  C.)  Grid  Points. 


VITA 


Frank  Broderick  Atkinson  vcas  born  on  13  \ovember 
1947  in  Nashville,  Tennessee  to  Robert  T.  and  Ruth  B. 
Atkinson.  After  graduation  from  Irving  Senior  High 
School  in  1965  he  attended  Arlington  State  College 
(now  the  University  of  Texas  at  Arlington) .  During 
his  studies  there  he  participated  in  the  cooperative 
education  program  as  a  Weight  Control  Engineer  at  Bell 
Helicopter  Company  in  Hurst,  Texas.  Following  gradua¬ 
tion  in  1970  with  a  Bachelor  of  Science  degree  in  Aero¬ 
space  Engineering,  he  enlisted  in  the  U.  S.  Air  Force  as 
a  Ground  Radio  Communications  Eguipment  Repairman.  He 
subsequently  attained  the  rank  of  Sergeant  whereupon  he 
was  selected  to  attend  Officer’s  Training  School  (OTS) 
at  Lackland  AFB,  Texas.  Upon  graduation  from  OTS  in 
January,  1975,  he  was  assigned  to  the  3246  Test  Wing, 

Guns  and  Fuzes  Division,  Eglin  AFB,  Florida  as  a  Wing 
Munitions  Test  Engineer.  During  his  tenure  at  Eglin 
AFB,  he  was  selected  three  times  as  his  division’s 
nominee  for  the  Directorate  of  Test  Engineering  Test 
Engineer  of  the  Quarter  Award.  He  also  represented  the 
Directorate  of  Test  Engineering  as  the  nominee  for  the 
Lt.  Robert  L.  Sullivan  Award  for  the  most  outstanding 
junior  officer.  Capt.  Atkinson  was  assigned  to  the 
School  of  Engineering  of  the  Air  Force  Institute  of 
Technology  in  June  1979  in  the  Graduate  Aeronautical  Engi¬ 
neering  Program. 
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